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INTRODUCTION 


This  report  describes  a  computational  method  for  calculating  supersonic 
lnvlscid  flow  within  an  inlet  using  an  extension  of  the  SWINT  code  which  is 
described  in  References  1,  2  and  3.  The  SUINT  code  is  designed  to  calculate 
3-D  external  supersonic  flow  fields  on  a  missile  type  configuration.  It 
marches  a  known  nosetip  flow  field  down  the  length  of  the  missile  body.  The 
computational  domain  is  bounded  on  the  Inside  by  the  body  surface  and  on  the 
outside  by  the  bow  shock  which  is  tracked.  It  is  restricted  to  flow  fields 
which  are  supersonic  everywhere  and  is  specifically  designed  to  treat  thin 
lifting  surfaces.  The  points  interior  to  the  computational  domain  are 
described  using  a  weak  conservation  form  of  the  Euler  equations  while  a 
characteristic  analysis  is  applied  to  determine  the  relations  applicable  at 
the  body  fin  or  shock.  The  MacCormack  explicit  method  is  used  to  advance  the 
flow  field. 

The  modifications  to  SWINT  described  in  this  report  primarily  consist  of 
replacing  the  shock  relations  with  solid  wall  conditions.  The  computational 
domain  is  now  bounded  along  the  lower  and  upper  edges  by  the  body  and  cowl 
respectively.  The  inlet  computation  is  started  using  the  flow  field  at  the 
inlet  face  determined  by  the  SWINT  code.  This  flow  field  is  re-grid  at  the 

e 

inlet  face  to  exclude  portions  of  the  flow  field  outside  the  cowl  using  the 
Interface  program  COWLI  which  is  also  described  in  this  report. 

The  extended  SWINT  code  is  only  applicable  to  the  supersonic  portion  of 
the  inlet  and  falls  when  the  axial  Mach  number  becomes  less  than  unity. 
Internal  shocks  are  not  tracked  but  instead  captured  by  the  numerical 
scheme.  The  modified  SWINT  code  is  most  applicable  to  cylindrical  inlets  (not 
necessarily  circular)  since  the  left  and  right  edges  of  the  computational 
domain  are  treated  using  either  a  symmetry  or  antisymmetry  condition 
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implemented  in  cylindrical  coordinates.  However,  reflection  boundary 
conditions  may  be  used  to  simulate  planar  walls  that  are  aligned  with  the  body 
axis.  This  makes  it  possible,  in  principle,  to  treat  certain  non-cyclindrical 
cases.  Planar  flows  can  also  be  approximated  by  using  large  body  and  cowl 
radii. 

In  the  remainder  of  this  report  the  extension  of  SWINT  to  handle  inlets 
is  outlined  and  several  test  cases  are  presented.  The  interface  program, 
COWLI,  is  also  described.  This  program  calculates  induced  inlet  drag, 
recovery  pressure  and  mass  capture  as  well  as  rezoning  the  inlet  face  flow 
field.  Appendices  A,  B  and  C  present  a  listing  of  COWLI,  an  update  deck 
listing  for  converting  the  original  SWINT  to  the  extended  version  described  in 
this  report,  and  user  instructions  for  applying  SWINT.  A  sample  inlet  run  is 
provided  in  Appendix  D. 
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MODIFICATIONS  TO  THE  SWINT  CODE 


COWL  BOUNDARY  CONDITIONS 

The  major  modification  of  SWINT  involves  implementation  of  surface 
boundary  conditions  at  the  outer  edge  of  the  computational  domain.  This 
development  parallels  that  for  the  body  surface  which  is  outlined  in 
References  1-3  and  given  in  detail  in  Reference  4.  For  completeness  the 
results  and  analysis  are  summarized  here.  The  physical  and  computational 
coordinates  applicable  to  an  inlet  configuration  are  shown  in  Figures  1  and  2 
respectively. 

On  the  cowl  surface  (X  ■  1)  the  normal  velocity  component  is  zero  which 
Implies : 

u  -  c  w  -  (c/c)  v  -  0  (1) 

z  <j> 

This  condition  is  supplemented  with  certain  characteristic  compatibility 

relations  associated  with  the  Euler  equations.  It  is  found  that  there  are 

three  Independent  characteristic  relations  which  are  admissible  on  X  -  1. 

These  can  be  written  as  a  system  of  quasi- linear  first  order  partial 

differential  equations  on  X  *  1  for  advancing  F  “  ln(p),V,  *  u(£/c)+  v  and  s. 

4  9 

The  resulting  relations  are 

H  ’  [XrX'  lx  ”  {pwt  X-  H  -  (a7w  +  a4v>l  +  P>]  ?  (2a) 


where  p  -p 1-  +  WA_ei  +  t  • 


(2b) 


+  pU2(Tg5Yz  +  %>  +  =  ^  (Tg^H-  c^/c)] 
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X  - 


a  (ercz) 
w2-a2  ‘ 


v  l  -  U  +(-^-)2) 


a7  ‘  TT  -  '  V  A 


a(c  /c)  c  c 

- $ - i[c  -  — L-±  -  (c 

az  ~  4 


-  S2/S  )Y  /Y  l 
4>  z  $ 


I  - 

5j  -  wX_(2-V2Kl/p)  ,  i2  -  cz-  X_  +  w2KjX_/p 

C3  -  WUK 1X_/p-l  ,  -  WVK1X_/p  +  Zfc 

*2  "  (|£)p  ■  ^"f^p  CiCj  -  -  p/h  for  a  perfect  gas) 

A 

Tg  “  *YY/gY  »  Tg  “  8zy/8y  *  G  "  GJ/r 
5  6 

v  -  1  +  (c  /c)  +  c 

a)  9  z 

A-Xu+|$v  +  Xw 
r  r  z 

A 

a^u  +  v/pw  (3a) 


•||  -  pvw  c  /c  -  p(T  Y  +  Z.M/c 
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nl  "  V4  *  n2  "  0  »  n3  “  »  n4  -  -1 

3s  B  3s 
3Z  "  '  v  3Y 


where  B  ■*  Yi^v  +  y  w 
r  z 

A 

Alternative  expressions  for  p  and  v  of  (2b)  and  (3b)  which  often  give 
improved  results  are: 


p  -  ?r  V.  +  [c  Y  +  1  c  Y  +  (Bw  -  YJ 

c  4  12  H  “7  z-3Y 


3P 


(2c) 


-  pB  (a5w  +  a^v) 


+  p-a.Kt,  S[t  ♦  ?]) 

*6  *6  S  f7  c  ♦  3Y  z 


where 


a5 


*v 

3c 

—  5  ■  c  /  Y 
3Y  zf  > 


3  (c  /c) 
*3  "  Ty 


~  (c^/c) 2)Y 
c  4> 


T  + 


fzx+  Vyx 


<bz  “  Cz> 
(c  -  b)  ' 


-  V**  |  (%  ~  V 

7  £X  (c  -  b) 


and 


v 


3V. 

<“  *  V  -  8T>  - 


1 

>+> 

c 


y  3P  pvw  c 
>  3Y  "  - 


(3c) 


11 


NSWC  TR  83-428 


Many  configurations  of  interest  have  sharp  corners  or  edges  such  as  those 
found  on  biconics  and  other  segmented  shapes.  If  the  upstream  cowl  surface 
velocity  normal  to  this  edge  is  supersonic,  either  a  shock  wave  or  an 
expansion  fan  will  be  attached  to  it  producing  a  discontinuity  in  the  surface 
flow  variables.  To  handle  this  situation,  an  oblique  shock  or  Prandtly-Meyer 
expansion  is  applied  at  the  edge  as  is  described  in  References  1  to  4.  In  the 
interior  these  discontinuities  are  captured  using  the  dissipative  and 
conservation  properties  of  the  interior  point  scheme.  Analogous  procedures 
are  applied  on  the  body  and  fin  surfaces. 
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COWL  INTERFACE  PROGRAM 

The  cowl  interface  program  operates  on  the  inlet  face  flow  field 
generated  by  an  external  SWINT  run.  It  is  designed  for  cylindrical,  but  not 
necessarily  circular,  inlets  and  rezones  the  flow  field  to  lie  within  the 
inlet.  In  addition,  this  program  calculates  recovery  pressure  throughout  the 
inlet  plane,  average  recovery  pressure  for  the  inlet  plane  and  the  flow 
entering  the  inlet,  mass  captured  by  the  inlet  and  induced  forces.  A  listing 
of  this  routine  is  provided  in  Appendix  A  and  user  instructions  are  outlined 
in  Appendix  C. 

The  average  inlet  plane  recovery  pressure  ratio,  average  inlet  rec<N~  'y 
pressure  ratio  and  the  mass  captured  by  the  inlet  are  determined  from 
Equations  (4a),  (4b)  and  (4c)  respectively: 


}{J«|  ‘  «l*p  J  /  p  t1  +  ^  "  ] 


y/(y-D 

rdrd$  (4a) 


2tt  c 


(  P  /p  )  •  — 

*t  *t®  Inlet  Pt 


J  H1 


fv-l) 

1  +  „2  rdr<4  (4b) 


J 


r  r 


pwrdrd$ 


(4c) 
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where: 


v  y  / v/  / 


rdrd<|i 


The  induced  forces  are  those  produced  by  the  action  of  pressure  on  the 
stream  surface  which  intersects  the  cowl  lip.  They  can  be  directly  determined 
by  storing  the  flow  field  upstream  of  the  inlet  face  and  tracking  the  stream 
surface  Intersecting  the  cowl  lip  back  through  the  flow  field  until  it 
intersects  the  bow  shock.  With  the  geometry  of  this  stream  surface  known,  the 
surface  pressure  along  it  can  be  integrated  to  produce  induced  drag  and 
lift.  For  complicated  bodies  at  incidence,  the  stream  surface  intersecting 
the  cowl  may  exhibit  a  very  complex  shape  and  this  type  of  procedure  is  both 
laborious  and  difficult  to  Implement.  An  alternative  approach  is  to  balance 
forces  and  moments  acting  in  the  control  volume  illustrated  in  Figure  3.  Here 
the  induced  forces  are  determined  by  performing  integrations  at  the  inlet  face 
plane.  The  resulting  equations  for  axial,  normal  and  yaw  force  are: 


■/./ 


p 

[pwCw^  -  w)  +  -r~-  -  p]  rdrd<j) 
^00  00 


/*  c 

J 


pw  [u  -  u]  rdrd$ 
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I 


I 


pw  [v^  -  v]  rdrd$ 


(5c) 


where  u  *  -  cos$u  +  sin$v 
v  -  sln$u  +  C08$V 

Figure  4  Indicates  the  sign  convention  for  these  forces.  Although  this 
procedure  Is  much  simpler  than  the  first  approach  described,  the  results  aust 
be  carefully  scrutinized.  There  may  be  cases  where  the  values  of  the 
Integrands • In  (5)  are  the  small  differences  between  two  large  nuabers.  Both 
approaches  have  been  tested  on  sample  axlsyaaetrlc  cases  and  results  agree 
within  several  percent.  An  additional  check  on  the  accuracy  of  this  procedure 
can  be  accomplished  by  comparing  the  force  coefficients  calculated  by  SKINT 
with  those  of  equation  (5).  Here  the  SKINT  calculated  force  coefficients 
correspond  to  those  acting  on  the  portion  of  the  body  forward  of  the  Inlet 
plane,  and  equations  (5)  are  Integrated  from  the  body  to  the  shock.  To  adjust 
for  the  assumption  used  in  SKINT  that  the  based  pressure  Is  pn  ,  A* A,  must 
be  added  to  the  drag  force  coefficient  calculated  by  SKINT.  Here  A*  is  the 
center  body  cross-sectional  area  at  the  Inlet  plane.  The  drag  force 
coefficients  computed  by  these  two  different  techniques  generally  agree  to 
within  2Z.  Discrepancies  between  the  SKINT  calculated  normal  and  yaw  force 
and  the  results  from  (Sb)  and  (5c)  are  greater,  with  errors  of  approximately 
I5Z  and  5Z  occuring  on  a  cone  at  5°  Incidence  at  Mach  2  and  Mach  4 
respectively.  In  order  to  provide  a  guide  to  the  accuracy  of  the  calculated 
induced  forces,  the  interface  program  COKLI  computes  this  comparison  for  each 
force  component. 
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RESULTS 


The  extended  SWINT  code  has  fc  :n  applied  to  a  number  of  different 
configurations,  three  of  which  are  Illustrated  In  this  section.  In  each  case, 
measured  surface  pressures  are  compared  to  calculated  values.  The 
computations  were  carried  out  using  19  points  between  the  body  and  the  cowl. 

In  the  axisymmetric  cases  three  circumferential  planes  were  used,  while  in  the 
other  situations  this  number  was  increased  to  13.  The  tested  inlets  were  of 
the  mixed  compression,  asymmetric  type  with  a  translating  center-body.  All 
cases  feature  boundary  layer  bleed  and/or  slots  to  reduce  the  thickness  of  the 
boundary  layer. 

The  first  inlet  considered  is  described  in  Reference  3  and  results,  along 
with  a  sketch  of  the  inlet,  are  shown  in  Figure  5.  This  example  features  a 
free-stream  Mach  number  of  2.3,  0°  incidence,  boundary  layer  bleed  and  a 
center  body  scoop  upstream  of  the  throat.  The  geometry  was  approximated  using 
a  piece-wise  continuous  function  generated  from  a  tabular  listing  of  the 
centerbody  and  cowl  profiles  which  was  provided  in  Reference  6.  Derivatives 
were  approximated  using  a  central  differencing  of  the  surface  locations 
evaluated  with  the  local  computation  step  size,  ±Az.  The  scoop  upstream  of 
the  throat  was  simulated  using  the  inlet  option  of  SWINT  which  excludes  from 
the  calculation  that  portion  of  the  flow  field  entering  the  scoop. 

Figure  6  illustrates  a  comparison  between  SWINT  results  and  those 
measured  in  Reference  6.  This  inlet,  which  is  depicted  in  Figure  6,  features 
wall  bleed,  a  free-stream  Mach  number  of  2.5,  and  an  incidence  of  C P  .  The 
body  and  cowl  geometries  are  described  using  the  cubic  splines  provided  in 
Reference  6. 
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The  final  case  considered  is  illustrated  in  Figure  7  and  features  an 
inlet  at  Mach  3.3  and  incidence  of  3°.  The  experimental  data  used  for 
comparison  was  generated  in  Reference  7  and  reported  in  References  7  and  8. 
Results  from  the  leeward  and  windward  planes  are  compared  with  experiment. 

The  geometry  description  was  generated  from  a  tabular  listing  of  the  body  and 
cowl  profiles.  Central  differencing  was  again  used  to  generate  needed  body 
and  cowl  derivatives. 

The  results  shown  in  Figures  6,  7  and  8  are  in  reasonable  agreement  with 
experiment.''  However,  the  predicted  location  at  which  shocks  strike  the 
centerbody  or  cowl  is  downstream  of  the  measured  one.  This  is  to  be  expected 
since  the  effective  distance  between  the  centerbody  and  the  cowl  is  decreased 
by  the  presence  of  the  boundary  layer. 

In  some  of  the  calculated  cases,  it  is  possible  to  march  through  the 
diffuser  to  the  end  of  the  inlet  without  encountering  subsonic  flow.  In  cases 
where  the  throat  Mach  number  is  greater  than  unity,  the  exit  conditions 
determine  the  location  of  the  terminal  shock  and,  hence,  subsonic  regions  in 
the  inlet.  The  marching  method  currently  employed  precludes  application  of 
the  downstream  boundary  conditions,  and  the  resulting  solution  represents  exit 
conditions  with  sufficiently  low  pressure  to  permit  supersonic  flow  throughout 


the  inlet 
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CONCLUDING  REMARKS 

The  SWINT  code  has  been  extended  to  allow  invlscid  calculations  to  be 
performed  on  the  supersonic  portions  of  inlets.  This  procedure  replaces  the 
bow  shock  tracking  procedure  with  solid  surface  boundary  conditions.  In 
addition,  an  interface  program  has  been  developed  which  rezones  the  external 
flow  field  upstream  of  the  inlet  face  to  Include  only  that  portion  of  the  flow 
field  entering  the  inlet.  The  Interface  program  also  calculates  forces,  mass 
captured  by  the  inlet  and  average  recovery  pressure  for  the  flow  entering  the 
inlet.  The  external  flow  field  upstream  of  the  inlet  can  be  determined  with 
either  the  original  version  of  SWINT  or  the  extended  version  described  in  this 
report.  The  extended  SWINT  code  is  best  suited  for  calculating  cylindrical 
inlets  (not  necessarily  circular)  and  is  restricted  to  geometries  where  the 
inlet  lip  lies  in  a  plane  perpendicular  to  the  missile  axis.  Comparisons 
between  calculation  and  experiment  have  been  performed  for  several 
axisymmetric,  external  compression  inlets  with  boundary  layer  bleed.  Computed 
surface  pressures  are  in  reasonable  agreement  with  experiment. 


COWL 


FIGURE  3.  CONTROL  VOLUME  FOR  CALCULATING  INDUCED  FORCES 
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COWL  STATIC  PRESSURE 

SST|  INLET,  MACH  -  2.30,  ANGLE  -  0.0 


BOW  SHOCK 
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FIGURE  5.  COMPARISON  OF  MEASURED  AND  CALCULATED  SURFACE  PRESSURES. 
EXPERIMENTAL  DATA  FROM  REFERENCE  5 
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CENTERBODY  STATIC  PRESSURE 
SST  INLET,  MACH  -  2.30,  ANGLE  -  0.0 


BOW  SHOCK 


FIGURE  6.  COMPARISON  OF  MEASURED  AND  CALCULATED  SURFACE  PRESSURES. 
EXPERIMENTAL  DATA  FROM  REFERENCE  6 
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FIGURE  6.  CONTINUED 
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COWL  STATIC  PRESSURE 

PRESLEY  INLET,  MACH  -  3.30,  ANGLE  -  3.0 
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FIGURE  7.  COMPARISON  OF  MEASURED  AND  CALCULATED  SURFACE  PRESSURES 
EXPERIMENTAL  DATA  FROM  REFERENCES  7  AND  8 
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CENTERBODY  STATIC  PRESSURE 
PRESLEY  INLET,  MACH  -  3.30,  ANGLE  -  3.0 
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FIGURE  7.  CONTINUED 
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speed  of  sound 

Location  of  the  body  surface 

Location  of  bow  shock 

Location  of  cowl  surface 

Enthalpy 

Stagnation  enthalpy 
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(r,4>,z) 


(u,v,w) 


In  (p) 


pressure 


recovery  pressure 


Velocity  vector 


Cylindrical  coordinates  (see  Figure  1) 


Entropy 


Velocity  components  in  cylindrical  coordinates  (see  Figure  1) 
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Computational  coordinates 


Computational  marching  step 


Density 
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a 

APPENDIX  A.  COWLI  PROGRAM  LISTING 

— 

1 

PROGRAM  COWL  1 NT  1 1 NPUT >64 • OUTPUT >6* • T APES* INPUT « T APE6*0UTPUT , 

COWL1NT 

2 

1  TAPE3«512.TAPE11) 

COWLINT 

3 

C. 

. THIS  PROGRAM  REZONES  THE  SHOCK  LAYER  AT  ThE  INLET  PLANE 

COWLINT 

4 

C 

TO  INCLUDE  ONLY  THE  F LOW  ENTERING  THE  INLET.  IT  CAN  ALSO 

COWLINT 

5 

* 

5 

c 

BE  USED  TO  6ENERATE  A  STARTING  FLOW  FOR  EXTERNAL  CALCULATIONS 

COWLINT 

6 

- 

c 

DOWN  STREAM  OF  THE  INLET  PLANE  WHEN  THE  BOH  SHOCK  OR  PORTIONS 

COWLINT 

7 

■ 

c 

OF  THE  BOW  SHOCK  ARE  INSIOE  THE  INLET. 

COWLINT 

B 

■ 

c 

INPUT  1 

COWLINT 

9 

c 

TAPE11  -  SWINT  RESTART  TAPE  AT  INLET  FACE 

COWLINT 

10 

10 

c 

BNEW. BZNEW, BPHNEW  -  NEW  BOOT  SHAPE  AT  INLET  FACE.  NEED  NOT 

COWLINT 

11 

c 

ALWAYS  BE  SPECIFIED. (SEE  1B00Y  OESCRIPTIO 

COWLINT 

12 

c 

CNEW.CZNEW.CPHNEW  -  SHAPE  OF  OUTER  BOUNDARY.  NEED  NOT  ALWAYS 

COWLINT 

13 

'  •  ' 

c 

BE  SPECIF IEO (SEE  ICOWL  DESCRIPTION! 

COWLINT 

1A 

* 

c 

IBODY  -  0  t  FLOW  AT  INNER  BOUNOARY  IS  NOT  ALTERED. 

COWLINT 

15 

15 

c 

BNEW, BZNEW .BPHNEW  NEED  NOT  BE  SPECIFIED. 

COWLINT 

16 

m 

c 

i  <  flow  at  inner  boundary  is  turneo  tangent  to 

COWLINT 

IT 

—  - 

a 

c 

THE  SURFACE.  BNEW .BZNEW, BPHNEW  MUST  BE  SPECIFIED. 

COWLINT 

IB 

c 

2  i  INNER  BODY  FLOW  QUANITIES  ARE  PRESCRIBED  ALONG 

COWLINT 

19 

c 

EACH  M  PLANE.  BNEW .BZNEW, BPHNEW  MUST  BE  SPECIFIED 

COWLINT 

20 

20 

c 

ICOWL  MUST  BE  EOUAL  0  AND  DOZ  MUST  BE  PRESCRIBED. 

COWLINT 

21 

c 

ICOWL  *  0  t  FLOW  AT  OUTER  BOUNOARY  IS  NOT  ALTERED. 

COWLINT 

22 

c 

CNEW.CZNEW.CPHNEW  MAY  BE  SPECIFIED.  OTHERWISE 

COWLINT 

23 

c 

OLO  VALUES  ARE  USED. 

COWLINT 

2A 

c 

1  >  FLOW  AT  OUTER  BOUNDARY  TURNED  TANGENT  TO  SURFACE 

COWLINT 

25 

25 

c 

CNEW.CZNEW  AND  CPHNEW  MUST  BE  SPECIFIED. 

COWLINT 

26 

■*> 

IL 

c 

2  !  OUTER  BOUNDARY  IS  A  MACH  CONE.  CNEW. CPHNEW  MUST 

COWLINT 

2T 

c 

BE  SPECIFIED.  CNEW  MUST  BE  GREATER  THEN  C  FOR  ALL 

COWLINT 

20 

T 

c 

PLANES. 

COWLINT 

29 

c 

AREA  -  REFERENCE  AREA!  OEFAULT  IS  BODY  CROSSECTIONAL  AREA 

COWLINT 

30 

30 

c 

at  the  inlet  entrance  plane. 

COWLINT 

31 

c 

RCLUST  -  R  DIRECTION  CLUSTERING!  DEFAULT  IS  UNIFORM  MESH. 

COWLINT 

32 

c 

DOZ  -  DISTANCE  FROM  COWL  LIP  TO  START1N6  PLANE.  ONLY  NEEDED 

COWLINT 

33 

c 

FOR  IB00Y*2. 

COWLINT 

3A 

c 

IPRINT  -  0  !  00  NOT  PRINT  FLOW  FIELD. 

COWLINT 

35 

» 

& 

35 

c 

1  S  PRINT  FINAL  FLOW  FIELD.  (OEFAULT) 

COWLINT 

36 

m 

c 

2  1  PRINT  ORIGINAL  AND  FINAL  FLOW  FICLO. 

COWLINT 

37 

w»  .  hJ 

c 

3  1  PRINT  ORIGINAL. FINAL  FLOW  FIELD  AND  JUMP  MESSAGE 

COWLINT 

30 

c 

OUTPUT! 

COWLINT 

39 

c 

TAPES  -  RESTART  TAPE  FOR  SWINT 

COWLINT 

40 

AO 

COMMON  NC.MC.K, IPRINT, PINF.DINF.PHIO. PI, RAD. Z.BZZ. GAMMA. M0T2.BMAX. 

CONST 

2 

1S1.S2.C1 ,C2, CONVR.PT INF  » 

CONST 

3 

1  C(100)  ,CZ(100)  .CPHI  (100)  .R  (25, 100). 0(25. 100)  >P(25.100)  «U(25.100)  . 

CONST 

4 

,  •„ 

1  V(25,100) ,W(25.100) .PHI (100) 

COMST 

5 

z 

2  . BNE W ( 1 0 0 > , CNE W (1 00 ) . DUMV ( 26 ) . DUMP ( 26 ) . OUHU ( 26 ) . DUMW ( 26 > . DUMD ( 26 > 

CONST 

6 

p 

■ 

♦5 

3  . BZNEW  (100)  .BPHNEW(IOO)  .CZNEW(IOO)  .CPhNEwUOO)  .R0L0(26) 

COMST 

T 

A  ,PTR(25.100> .RCLUST (100). PHIL (102) 

COMST 

0 

COMMON/RGASS/AX.HX.GX 

RGASS 

2 

REAL  mx.my.mz.mxz.myz.mzz 

COWLINT 

A3 

NAMEL IST/INPUTS/BNEW, BZNEW, BPHNEW. CNEW.CZNEW. CPHNEW 

COWLINT 

44 

50 

1  ,I800V, ICOWL, AREF. RCLUST, ODZ, IPRINT 

COWLINT 

A5 

c. 

....INITIALIZE  OATA 

COWLINT 

46 

DO  5  H*1 ,25 

COWLINT 

AT 

PNEw(M)*-l. 

COWLINT 

AO 

CNEW  (M )  *- 1 . 

COWLINT 

49 

55 

5  CONTINUE 

COWLINT 

50 

002*0. 

COWLINT 

51 

IPRINT  *  1 

COWLINT 

52 

PTRTOT  ■  0.0 

COWLINT 

53 

PTRINL  *  0.0 

COWLINT 

5A 

60 

AINL  *  0.0 

COWLINT 

55 

* 

AREAT  *0.0 

COWLINT 

56 

FAOOA  ■  0.0 

COWLINT 

57 

A 

FAOON  *  0.0 

COWLINT 

50 

K_ 

FADDY  *  0.0 

COWLINT 

59 

65 

FAOOAT*0. 

COWLINT 

60 

FAOONTsQ. 

COWLINT 

61 

- 

FADDYT*0. 

COWLINT 

62 

SUMINMaO. 

COWLINT 

63 

ICOWL  ■  1 

COWLINT 

6A 

TO 

IBODY*0 

COWLINT 

65 

.‘.V 

AREF  «  0.0 

COWLINT 

66 

c 

COWLINT 

67 

SL 

c. 

. ...REAO  FLOW  FICLO 

COWLINT 

60 

— 

c 

COWLINT 

69 

75 

R£AD( 1 l ) NC.MC, ATTACK. YAW, ACM.GAMM A, PINF.OINF, PHI O.K,Z, 

COWLINT 

TO 

2  FN.FY , FA, MX, MY ,MZ,FNZ,FYZ,FAZ,MXZ,NYZ,MZZ. 

COWLINT 

71 

.•  • 

3  (PHI (M) ,C(M> ,CZ (M) ,CPHI (M) ,M*1,MC> 

COWLINT 

72 

A  ,  ( (R(N.M) ,U(N,M> ,V(N,M) .WIN.M) ,P(N.M) .O(N.M) • 

COWLINT 

T3 

»- 

A-l 

-- 

- ■  ■  v  ■  " .  j 

~4 

1 
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5  M  ■  1 »HC) ,N*1 .NC) 

cowlint 

T4 

80 

irtEOFU  11)1000. 25 

COWLINT 

75 

.  ’  J 

c 

COWLINT 

76 

^ 

E 

25  CONTINUE 

cowlint 

77 

■ 

00  3  N*1»NC 

COWLINT 

78 

RCLU5T (N) * (N-l . ) / (NC-1 • ) 

COWLINT 

79 

85 

3  CONTINUE 

cowlint 

80 

c 

COWLINT 

81 

C . READ  GEOMETRY  DATA 

COWLINT 

82 

»■ 

READ (5. INPUTS) 

COWLINT 

83 

v  ■* 

► 

WRITE (6. INPUTS) 

COWLINT 

84 

90 

PI  «  4.0  •  ATAN (1*0) 

COWLINT 

85 

•1  ”  ; 

| 

CONVR*180./PI 

COWLINT 

86 

1 

P 

ATWO«0. 

COWLINT 

87 

MP?«MC«2 

COWLINT 

88 

MP1«M01 

cowlint 

89 

95 

00  10  M-2.MP1 

COWLINT 

90 

PHIL(N)*PHI(N-1> 

COWLINT 

91 

10  CONTINUE 

COWLINT 

92 

PHILID-PMI  ID 

COWLINT 

93 

PHIL (MP2) »PMl (MC) 

COWLINT 

94 

p  100 

IFIABS(PH10-2«PI> .GT.1.E-061GO  TO  11 

COWLINT 

95 

■g 

PHIL(1)*-PH1<2> 

COWLINT 

96 

PHIL! HP2 ) *2  *  *P 1 

COWLINT 

97 

11  CONTINUE 

COWLINT 

98 

4 

I 

. . COMPUTE  REFERENCE  AKEA 

COWLINT 

99 

[  105 

00  *  M  ■  1.  MC 

COWLINT 

100 

IM1  «  0 

COWLINT 

101 

1  ' 

IMC  ■  2 

COWLINT 

102 

■ 

ATWO  m  ATNO  ♦  0.25  •  (PHIL (M« IMC) -PHIL (M-IM1) 1  •  R 1 1 .M) -R 1 1 ,M> 

COWLINT 

103 

L. 

*  CONTINUE 

COWLINT 

104 

_  - * 

■  no 

ATNO  *  ATWO  •  2.0  •  PI  /  PHIO 

COWLINT 

105 

l  . 

m 

2  CONTINUE 

COWLINT 

106 

IF ( AREF .EO.O.O) AREF* AT WO 

cowlint 

107 

C . CALCULATE  CONSTANTS 

cowlint 

108 

6X*GAMMA 

COWLINT 

109 

115 

NCM1-NC-1 

COWLINT 

110 

RAD-PI/180. 

COWLINT 

111 

ALPHA* ATTACK-RAD 

COWLINT 

112 

S**SIN(ALPHA) 

COWLINT 

113 

vinf*sqrtigamma*pinf/oinf)«ach 

COWLINT 

114 

1 

P  12° 

YAWR  *  YAW  •  RAO 

COWLINT 

115 

r 

WINF  *  VINF  •  COS (YAWR)  •  COSIALPHA) 

COWLINT 

116 

1 

VI  INF  -  VINF  •  COS (YAWR)  •  S* 

COWLINT 

117 

1 

V2INF  «  -VINF  •  SIN(YAWR) 

COWLINT 

118 

-  ■-  : 

I 

PTINF  «  PINF*  (1.0  ♦  (GAMMA-1.0)  /  2.0  •  ACH—2.0)  ••  (GAMMA/ 

COWLINT 

119 

.■  .  • 

|  125 

1  (GAMMA-1.0)) 

COWLINT 

120 

.  . 

I 

HO*GAMMA*P INF/ ( 1  GAMMA-1 . ) *0 INF ) ♦ V I NF**2/2 . 

COWLINT 

121 

. 

t* 

CALL  RGAS(PINF.OINF.SINF) 

COWL INT 

122 

. 

C . PRINT  FREE  STREAM  CONDITIONS 

COWLIJDT 

123 

WR ITE (6. 5000 ) ACH. ATTACK. YAW. VINF. PINF tDINF. HO tSINF, PTINF 

COWLINT 

124 

—  -■  -M 

130 

5000  FORMAT ( 1H1 ,5X» 38H*****»*  FREE  STREAM  CONDITIONS  •••••**,/ 

COWLINT 

125 

1 .1SX.15HMACH  NUMBER  .F15.7./ 

COWLINT 

126 

2. 15X, 15HANGLE  OF  ATTACK ,F 15.7,/ 

COWLINT 

127 

3.15X.15HYAW  ANGLE  .F15.7./ 

COWLINT 

128 

4. 15X. 15HVINF  .F15.7,/ 

COWLINT 

129 

135 

5.15X.15HPINF  ,F 15.7 ,/ 

COWLINT 

130 

- 

6. 15X . 15H01NF  .FIS. 7,/ 

COWLINT 

131 

- 

7.15X.15HH0  .FIS. 7,/ 

COWLINT 

132 

| 

8* 15X. 15HSINF  .F15.7 ./ 

COWLINT 

133 

, 

9.15X.15HPTINF  .F15.7,///) 

COWLINT 

13* 

-  *i  “I 

1*0 

C . PROBLEM  SET  UP 

COWLINT 

135 

WRITE(6*5001)NC*MC.IBODY,DDZ«ICOVL 

COWLINT 

136 

WRITE (6.5003) (N.RCLUST(N) .N*1.NC> 

COWLINT 

137 

COWLINT 

138 

1  6X . 1 HN. 7X . 1 OHC LUSTER 1NG . / 

COWLINT 

139 

1*5 

2  (2X.15.F15.7)) 

COWLINT 

140 

WRITE (6.5005) 

COWLINT 

141 

- 

5005  FORMAT (///) 

COWLINT 

1*2 

L 

COWLINT 

1*3 

I  .15X.15HNC  .15./ 

COWL INT 

14* 

150 

2  .15X.15HMC  .15./ 

COWLINT 

1*5 

3  , 15X . 15HIR00Y  .15./ 

COWLINT 

1*6 

3  » 15X, 15H0DZ  ,F10. 5,/ 

COWLIN' 

1*7 

4  » 15X. 15HIC0WL  .15,///) 

COWLINT 

148 

H0T2-2.-H0 

COWLINT 

1*9 

155 

El  *  PINF  /  DINF  /  WINF 

COWLINT 

150 

E3  «  0.5  •  DINF  •  VINF-VINF  •  AREF 

COWLINT 

151 

C . MAKE  NEW  GEOMETRY  AXISYMETRIC  IF  ONLY  ONE  VALUE  IS  GIVEN 

COWLINT 

152 

DO  6  M*2.MC 

COWLINT 

153 

IF  (CNEW(M) .GT.O , ) GO  TO  7 

COWLINT 

15* 

100 

CNEW(M)«CNEW(M-1) 

COWLINT 

155 

CZNEW (M) *CZN£W (M-l ) 

COWLINT 

156 

A- 2 

> 

f  *  -  • 

i 
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• 

*' 

* 

CPMNEV (M) aCPHNEW CN-1 1 

COWLINT 

157 

’ 

7  CONTINUE 

COWL1NT 

158 

— - - 

IF(BNEW(M).GT.O.)GO  TO  6 

COWLINT 

159 

• 

165 

BNEW (M) «BNEW (M-l ) 

COWLINT 

160 

BZNEW (M) *BZNEW (M-l ) 

COWLINT 

161 

•  • 

BRMNEW IN) aBPHNEW (M-l 1 

COWLINT 

162 

•  ■  : 

6  CONTINUE 

COWLINT 

163 

C. 

. use  tape  read  geometry  if  no  value  is  given 

COWLINT 

16* 

.  •  '■ 

1T0 

00  8  Mal.MC 

COWLINT 

165 

• 

IF (BNEW (M) «LT .0.) BNEW (M) aRl 1 . M) 

COWLINT 

166 

•  % 

IF (CNEW(M) .GT. 0.160  TO  8 

COWLINT 

167 

.  .  ^ 

CN£W(M)*C(N) 

COWLINT 

166 

A 

CZNEW (M) *CZ (M) 

COWLINT 

169 

175 

CPHNEW (M) aCPHI (M) 

COWLINT 

170 

8  CONTINUE 

COWLINT 

171 

WRITE (6.5002) (M.BNEW(M) .BZNEW (Ml . BPHNEw (Ml ,CNEW(M) . 

COWLINT 

172 

■ 

1  CZNEW (M) .CPNNEW (N).NaJ.MC) 

COWLINT 

173 

.  ■ 

5002  FORMAT (*0**5X.29M*******  COOL  GEOMETRY  •••••••,/ 

COWLINT 

17* 

.  ■ 

ISO 

1  5A.1HM.12X.1NB.13X.2MBZ.1 1X.4HBPHI .1*X. 1HC. 13X.2HCZ. 1 1X.4HCPHI/ 

COWLINT 

175 

-T 

2  (•  •.IS.6F15.7I 1 

COWLINT 

176 

C. 

....PRINT  ORIGINAL  FLOW  FIELD 

COWLINT 

177 

IF (IPRINT.GT.l)WRITE (6.5006) 

COWLINT 

178 

• 

5006  FORMAT (1H1.10X. 19H0RI6INAL  FLOW  FIELD) 

COWLINT 

179 

165 

IF ( I PRINT ,6T . 1 ) CALL  OUTPUT (ACN. ATTACK. YAW) 

COWLINT 

160 

C 

COWLINT 

161 

COWLINT 

162 

c 

COWLINT 

183 

C. 

....PRE-COWL  FLOW  FIELD.  SMOCK  IS  OUTER  BOUNORY 

COWLINT 

184 

190 

c 

COWLINT 

185 

00  31  N  .  1,  AC 

COWLINT 

186 

IMO  >  1 

COWLINT 

187 

m 

INC  ■  2 

COWLINT 

188 

IFLA6  «  0 

COWLINT 

169 

195 

C . LOCAL  VALUES  ARE  CALCULATED 

COWLINT 

190 

C. 

. SOME  ARE  AREA  WEIGHTED  AND  SUMMED 

COWLINT 

191 

'  ■  ‘  * 

DO  30  N  a  1,  NC 

COWLINT 

192 

... 

AX  a  SORT I6X*P (N.M) /D (N.M) ) 

COWLINT 

193 

AHACM  a  SORT (U (N.M) **2.0  ♦  V(N.M)a»2.0  •  W (N.M) ••2.01 /AX 

COWLINT 

19* 

200 

PTR(N.N)  a  P (N.M)  •  (1.0  *  (OAMMA-1. 01/2.0  • 

COWLINT 

195 

1  AMACN**2,0)  ••  (GAMMA/1 GAMMA- 1,0)) 

COWLINT 

196 

m 

INI  a  1 

COWLINT 

197 

. 

INC  a  1 

COWLINT 

198 

IF  (N  .£6.  1)  INI  a  0 

COWLINT 

199 

205 

IF  (N  .EG.  NC)  INC  a  0 

COWLINT 

200 

’  •  . 

c. 

COWLINT 

201 

.■  s 

c. 

COWLINT 

202 

/  ■  * 

AREA  a  0.25  •  (PHIL (M*IMC> -PHIL (M-IMl ) )  •  ( (0.5* (R (N.M) ♦ 

COWLINT 

203 

1  R(N*INC.M) ) )**2.0  -  <0.5«fR(N-lNl»N)*R(N,N)))*»2.0) 

COWLINT 

20* 

210 

AREAT  a  AREAT  •  AREA 

COWLINT 

205 

E2  a  0  (N.M)  •  W  (N.M)  •  AREA 

COWLINT 

206 

9 — , 

VI  a  — U  (N.M)  •  COS  (PHI  (M)  )  «  VIN.M)  •  SIN(PHKM)) 

COWLINT 

207 

V2  a  U (N.M)  •  SIN (PHI (M) )  .  VIN.M)  •  COS(PHIIM)) 

COWLINT 

208 

FAOOAT  •  FAOOAT  ♦  ( (WIMF-W (N.M) )  »  El  -  P(M.M)/0(M.M) 

COWLINT 

209 

’  ■  4 

215 

1  /W (N.M) )  •  F2 

COWLINT 

210 

FAODNT  a  FADONT  .  (V1INF-V1)  •  E2 

COWLINT 

211 

FAOOYT  a  FADDYT  ♦  (V2INF-V2)  •  12 

COWLINT 

212 

PTRTOT  a  PTRTOT  .  PTR(N.M)  •  AREA 

COWLINT 

213 

, 

IF (R (N.M) .LE.CNEW (M) ) 60  TO  29 

COWLINT 

21* 

220 

C . ELEMENT  INTERSECTED  BY  COWL 

COWLINT 

215 

• 

IF  (N  .EQ.  1)  60  TO  29 

COWLINT 

216 

C. 

COWLINT 

217 

C. 

COWLINT 

218 

* 

IF  1 IFLAG  .to.  1)  GO  TO  28 

COWLINT 

219 

225 

C. 

COWLINT 

220 

C. 

COWLINT 

221 

AREAB  a  0.25  •  (PHIL (M. IMCI -PHIL (M-IMl ) 1  •  ( (0.5»(RIN,M) ♦ 

COWLINT 

222 

l  R(N-l.M) ) )**2«0  -  CNEW (M) ••2.0) 

COWLINT 

223 

NB  m  H 

COWLINT 

22* 

230 

IF  (0.5*(H(N.MI*R(N-1.M>>  .GT.  CNEW(MI)  NB  ■  N  -  l 

COWLINT 

225 

— « 

£2  a  O(NB.M)  •  W(NB.M)  •  AREAS 

COWLINT 

226 

VI  a  -U(NB.M)  •  COS  (PHI  (M) )  .  V(NB.M)  •  SINIPHKM)) 

COWLINT 

227 

V2  a  U(NB.M)  •  SIN (PHI (H) )  •  V(NB.M)  •  C0S(PHI (Ml ) 

COWLINT 

228 

* 

FADOA  a  FADDA  ♦  UWINF-W  (NB .Ml )  ♦  El  -  P (NB .M)  /D (NB.MI 

COWLINT 

229 

235 

1  /V(NR.M))  •  E2 

COWL INT 

230 

FAODN  a  FADON  «  (V11NF-V1)  •  E2 

COWLINT 

231 

FADOY  a  FADDY  «  (V2INF-V2)  •  E2 

COWLINT 

232 

SUM INM  a  SUM1NH  -  D(NB.M)  •  W (NB.MI  •  AREAB 

COWLINT 

233 

PTRINL  ■  PTRINL  •  PTR(NB.M)  •  AREAB 

COWLINT 

23* 

ft _ 4 

2*0 

AINL  a  AINL  -  AREAS 

COWLINT 

235 

1 

IFLAG  a  1 

COWLINT 

236 

.  .  .1 

C. 

COWLINT 

237 

-  ■ 

c. 

. OUTSIDE  COWL 

COWLINT 

236 

28  CONTINUE 

COWLINT 

239 

A-3 

1  _  H 
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2*5 

E2  *  0(N,M)  •  W (N.M)  •  AREA 

COWL  I NT 

2*0 

VI  *  -U(N.M)  •  COSIPHI(M))  «  V(N.M)  •  SIN(PHKM)) 

COWLINT 

2*1 

V2  ■  U(N.M)  •  SIN (PHI (M) )  ♦  V(NtM)  •  C0S(PHl(M>) 

COWLINT 

2*2 

FAOOA  ■  FAOOA  *  < (WINF-W(N.M) 1  ♦  El  -  P(N.N)/D(N, 

Ml 

COWLINT 

2*3 

1  /WIN. Ml)  •  E2 

COWLINT 

2** 

250 

FAOON  a  FADON  *  (V11NF-V1)  •  E2 

COWLINT 

2*5 

FADOY  a  FADDY  •  (V2INF-V2)  •  E2 

COWLINT 

2*6 

60  TO  30 

COWLINT 

2*7 

C 

. INSIDE  COWL 

COWLINT 

2*8 

29  SUNINH  ■  SUM INN  «  D (N.M)  •  W(N.M)  •  AREA 

COWLINT 

2*9 

255 

PTRINL  *  PTR1NL  «  PTR(N.M)  •  AREA 

COWLINT 

250 

AINL  *  AINL  •  AREA 

COWLINT 

251 

30  CONTINUE 

CONJ.INT 

252 

IF(R(NC.M).6T.CNEW(H))G0  TO  31 

COM.  I  NT 

253 

c 

. SHOCK  INSIDE  COWL 

COWLINT 

25* 

260 

AREA-.25* (PHIL (M« INC) -PHIL (M-IM1 ) >*(CNEH(N)**2-R(NC< 

HI *a2) 

COWLINT 

255 

suminm»suminm.dinf»vinf»area 

COWLINT 

256 

PTRINL *PTRINL»PTINF»AREA 

COWLINT 

257 

AINL* AINL 'AREA 

COWLINT 

258 

31  CONTINUE 

COWLINT 

259 

265 

c 

. TOTAL  FORCES 

COWLINT 

260 

FADOAT*FADDAT#2 . «P I /PM  1 0 

COWLINT 

261 

F AODNTaF  ADDNT *2 • #P I /PH 1 0 

COWLINT 

262 

IF(ABS(PMIO-2.»PI) .GT.1.E-06)FADDYT»0. 

COWLINT 

263 

IF(AeS<<PHIO-PI)*<PHIO-2.«PI)).0T.l.E-O6>FADDNTa0. 

COWLINT 

264 

2T0 

faa*fa*atvo*pinf 

COWLINT 

265 

ERRA* (FADOAT-FAA) /FAA*100 . 

COWLINT 

266 

ERRNa (FADONT-FN) /SIGN < AMAX1 (ABS(FN) . l.E~08) »FN)*100. 

COWLINT 

267 

ERRY«(FAOOYT-FY)/SIGN(AMA*l (ABS(FY) . l.E-OB) .FT) »100. 

COWLINT 

268 

c, 

. INDUCED  DRAG 

COWL INT 

269 

275 

FAD0A«FA0DA*2.*PI/ (PMIO*E3) 

COWLINT 

270 

FAODN«FADON*2.*PI/ (PMI0*E3) 

COWLINT 

271 

FAOOY«FAOOY/E3 

COWLINT 

272 

IF(ABS(PHI0-2.»PI) ,GT.1.E-06)FAODY*0. 

COWLINT 

273 

IFIABSI (PMI0-PI)»(PMI0-2.«PI) ) .GT.1.E-06IFADON-0. 

COWLINT 

27* 

280 

c, 

. TOTAL  PRESSURE  RECOVERY  OATA  PRINTED 

COWLINT 

275 

PR  a  PTRTOT  /  AREAT 

COWLINT 

276 

PRINL  «  PTRINL  /AINL 

COWLINT 

277 

AREAT>AREAT*2.*PI/PHI0 

COWLINT 

278 

AINL«AINL*2.*PI/PH10 

COWLINT 

279 

285 

SUM INM*SUMINM*2  ••PI/PMIO 

COWLINT 

280 

WHITE  (6.2000) 

COWLINT 

281 

2000  FORMAT (3AH1 INLET  PLANE  FLOW  FIELD  PARAMETERS) 

COWLINT 

282 

PRINLHaPRINL/PTINF 

COWLINT 

283 

PRRsPR/PT INF 

COWLINT 

28* 

290 

WRITE (6. 20101PRR .PRINLR. AREAT. AINL. SUHINM.FADDA.ERRA. 

COWLINT 

285 

1  FADDN.ERRN.FAODY.ERRY.AREF 

COWLINT 

286 

2010  FORMAT (1H0.10X.A4HSHOCK  LAYER  AVCRA6E  PRESSURE  RECOVERY 

RATIO  . 

COWLINT 

287 

1  FIS. 7./. 

COWLINT 

288 

1  1 OX. ASH  INLET  AVERAGE  PRESSURE  RECOVERY  RATIO 

•Fi5#7*/t 

COWLINT 

289 

295 

2  10X.45H  SHOCK  LAYER  CROSSECTIONAL  AREA 

.F15.T./. 

COWLINT 

290 

3  10X.A5H  INLET  ENTRANCE  CROSSECTIONAL  AREA 

•FIS. 7./. 

COWLINT 

291 

A  10X.A5H  MASS  CAPTURED  BY  THE  INLET 

COWLINT 

292 

*  10X.A5H  ADDAT1VE  AXIAL  FORCE  COEFFICIENT 

•F15.7, 

COWLINT 

293 

*  AX.29H  TOTAL  0RA6  ERROR  .F10.A.AH  0/0./. 

COWLINT 

29* 

300 

5  10X.A5H  AQOATIVE  NORMAL  FORCE  COEFFICIENT 

.F 15*7. 

COWLINT 

295 

5  4X.29M  TOTAL  NORMAL  FORCE  ERROR  .F10.4.AH  0/0./. 

COWLINT 

296 

6  10X.A5H  AOOATIVE  YAW  FORCE  COEFFICIENT 

•F15.7. 

COWLINT 

297 

6  AX.29H  TOTAL  YAW  FORCE  ERROR  .F10.A.AH  0/0./. 

COWLINT 

298 

7  10X.45H  REFERENCE  AREA 

♦FI5.7/1H1 > 

COWLINT 

299 

305 

DO  *5  Mal.MC 

COWLINT 

300 

COWLINT 

301 

ROLO (MC*1)*2.*R(NC.M)«R(NCM1.M) 

COWLINT 

302 

S1»SIN(PHI(M)> 

COWLINT 

303 

C1*C0S(PMI (M) I 

COWLINT 

30* 

310 

DUMP (NC*1 ) aPINF 

COWLINT 

305 

0UM0(NC«1)*0INF 

COWLINT 

306 

DUMU(NC*l)»-VINF»S*»Cl 

COWLINT 

307 

DOMV INC* 1 ) *VINF*SA*S1 

COWLINT 

308 

OUNWlNC*!) "SORT ( V1NF»VINF-DUMU (NC» 1 ) ••2-DUMV (NC» 1)»«2) 

COWLINT 

309 

315 

DO  50  Nal.NC 

COWLINT 

310 

C. 

COWLINT 

311 

ROLO(N) «R(N.N) 

COWLINT 

312 

DUMP  IN)  aP  IN  .Ml 

COWLINT 

313 

OUMO(N)aOIN.N) 

COWLINT 

31* 

320 

OUMU(N)autN.M) 

COWLINT 

315 

DUMV (Ml aV IN, Ml 

COWLINT 

316 

DUMW (N) aw (N.M) 

COWLINT 

317 

C, 

COWLINT 

318 

R(N.M)  a  SNEW(M)  .  (CNEW (M) -BNEw (Ml  1  a  RCLUST(N) 

COWL INT 

319 

325 

50  CONTINUE 

COWLINT 

320 

A- 4 


Jk. 
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00  60  N«1*NC 
00  65  J*1  .NC 

C.. ......... ...DETERMINE  IF  R(N,M|  IS  INSIDE  010  SHOCK 

IF< (R(N.N)-ROLO(JI I • (R (N.N) -ROLD ( J»1 > > .6T  . 0.160  TO  65 
JJ«J 

GO  TO  70 

65  CONTINUE 

. . OUTSIOE  OLO  SHOCK.  FREE  STREAM  <F>  VALUES  GIVEN  TO  RtN.N) 

P(N.M)aPINF 

D(N,M)bDINF 

U(N.M»«-VINF»S*»C1 

V(N.M)bV1NF»S4«S1 

V(N,H)bS0RT(VINFbVINF-U(N.M)*42-V(N,M)«#2I 
GO  TO  60 
70  CONTINUE 

. . INS10E  OLD  SHOCK.  INTERPOLATED  (I)  VALUES  GIVEN  TO  RtN.N) 

JJ1 *JJ* 1 

FAC«(R(N.M)-ROLO(JJ))/<ROLDtJJll-ROLDlJ>i>  > 

U(N.H)bOUMU(JJ)  ♦(OUHU(JJl)-DUMU(JJ)  MFAC 
V(N.M)bOUMV(JJ) . (OUNV(JJl)-OUMV(JJI IbFAC 
P (N.H) sOUMPI JJ) ♦ tOUHP  t JJl I -DUMP ( JJ>  >  »FAC 
0  <N.M) aOUHOt JJ) ♦ (OUHO ( JJI ) -DUHDt JJ) ) »FAC 
V (N.H) bSQRT (M0T2-P (N.H) »GAHHA*2./ t IGAHHA-l.) •OtN.N) ) 

1  -U(N,M)b*2-V(N»M> **2) 

60  CONTINUE 

IF ( IBODY.EQ.O)GO  TO  A1 

. . BOOT  JUMP  (VALID  ONLY  FOR  AN  UNSVEPT  LEADING  EDGE) 

BZ0*(U(1.M)-V(1*M) •BPHNEV(H) /BNEV (M))/V(1»H) 

OELBZbBZO-BZNEV (N) 

SI0E-ASINI1./ACH) 

IFIABS(DELBZ) .LE.1.E-06IG0  TO  A3 
OELBPbO. 

SIOE  ■  1. 

CALL  JUMPSTtOELBP.OELBZ.M.l.BNEVIH) .BZNEV (N) .BPHNEV (H) .SIDE) 
U(l.H)  ■  SIGN  (U(l.M).  BZNEV (M) I 
IF ( I BOOT .EG. 1 ) 60  TO  41 

. . UNIFORM  FLOM  FROM  BODY  TO  SHOCK 

43  CONTINUE 

CZNEH (M). TAN (SIOE.ATAN (BZNEV (M) )) 

CNEV (M) >BNEV (M) .DOZ» ( CZNEV (M> -BZNEV ( M) ) 

CPHNEV (M) aBPHNEV (H) 

Z«Z*00Z 

ICOVL-O 

00  42  N»2,NC 

R (N.H) >BNEV (M) « (CNEV (M) -BNEV (M) >  VRCLUST (N) 

U(N,HI»U(1.M) 

V(N.H)«V(1.H) 

V(N.H)>V(1.M) 

P(N,H)«P(l,M) 

0(N.H)-0(I.H) 

42  CONTINUE 

41  CONTINUE 

IF  (ICOML  .EO.  1)  GO  TO  44 

. . COMPUTE  CZ  FOR  MACH  ZONE 

IF ( ICOVL.CQ.2) 

1  CZNEV(H)»(U(NC.M) -CPHNEV (M)»V(NC.H) /CNEV (M) 

2  .SORT ( 6AHHA*P (NC.M) /O (NC.HI ) ) /V (NC.M) 

GO  TO  45 

44  CONTINUE 

C . .  JUMPIVALID  ONLY  FOR  UNSVEPT  LEADING  EDGE) 

BZOa(UtNC.M) -V (NC.M) «CPHNEV (H) /CNEV (M) ) /V (NC.M) 
DELBZaBZO-CZNEV(M) 

IF ( ABS(OELBZ) .LE. 1 .E-6) 60  TO  45 
DELBP-O. 

SIOE—1. 

CALL  JUMPST (OELBP.DELBZ.M.NC. CNEV (M). CZNEH IN). CPHNEV (N>. SIDE) 
U(NC.H)  b  SIGN  (U(NC.M).  CZNEV(H)) 

45  CONTINUE 

VRITE(3)NC.HC.ATTACK.YAV.ACHiGAMMA.PINF.DINF.PHI0.K.Z 

1  .fn.fy.fa.ha.my.hz.fnZ.fyz.faz.hxz.htz.hzz 

2  . (PHI (M) .CNEV (HI .CZNEV(M) .CPHNEV (M) .Mb  I, MCI 

3  . ( (R (N.H) .U(N.M) ,V (N.H) «V (N.H) .P(N.M) .O(N.M) .Mb).MC) .Mb), NC> 
C.... .PRINT  FINAL  FLOV  FIELO  DATA 

VRITE (6.5007) 

5007  FORMAT ( 1HI t 10K. 16HFINAL  FLOV  FIELD) 

IF (IPRINT.GT.Q)CALL  OUTPUT (ACH. ATTACK, YAH) 

STOP"  PROGRAM  START" 

1000  CONTINUE 

STOP"  TAPE  REA01N  ERROR" 

END 


COVLINT 

COHLlNT 

COVLINT 

COVLINT 

COVLINT 

COHLlNT 

COVLINT 

COHLlNT 

COVLINT 

COVLINT 

COHLlNT 

COVLINT 

COVLINT 

COHLlNT 

COHLlNT 

COVLINT 

COVLINT 

COHLlNT 

COHLlNT 

COHLlNT 

COHLlNT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COHLlNT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COHLlNT 

COVLINT 

COVLINT 

COHLlNT 

COHUNT 

COHLlNT 

COVLINT 

COHLlNT 

COML INT 

COVLINT 

COHLlNT 

COHLlNT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COVLINT 

COHLlNT 

COVLINT 

COVL INT 

COVLINT 

COVLINT 

COHLlNT 

COVLINT 

COVLINT 

COVL INT 

COVL INT 

COVLINT 

COVLINT 
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1 

SUBROUTINE  OUTPUT (ACH, ATTACK, YAW) 

OUTPUT 

2 

COMMON  NC*MC,K,IPRlNT,PINF,DINF,PHIO,PI , R AO, Z, BZZ, 6 AMMA, HO T2.BMAX, 

CONST 

2 

1S1.S2. Cl, C2.CONVR.PT INF, 

CONST 

3 

1  CUOOI.CZ  (100)  ,CPH1  (100)  ,R(2S,100)  .0(25,1 00)  ,P(25,100).  U  (25,100  >, 

CONST 

4 

s 

1  V(25. 100), 0(25,100). PHK100) 

CONST 

5 

2  ,BNEW (100) ,CNEW (100) , DUMV (26) «OUMP (26) . DUMU (26) , DUMB ( 26 ) , OUMD ( 26 ) 

CONST 

6 

3  .BZNEW (100) , BPHNEW (100) «CZNEW (100) »CPhNEW (100) ,ROLO(26) 

CONST 

7 

4  , PTR(25*100> .RCLUST (100) .PHIL (1 02 » 

CONST 

B 

C0MM0N/R6ASS/AX , HX ,6X 

RGASS 

2 

10 

WRITE  (6,3000)  ACM, ATTACK, YAW 

OUTPUT 

5 

3000  FORMAT (•  MACH  NO  IS*,1PE1S.7,5X,*ANGLE  OF  ATTACK  IS*,1PE1S.7,SX, 

OUTPUT 

6 

1  •  ANGLE  OF  SIOESLIP  IS*,1PE15.7> 

OUTPUT 

7 

00  100  M>1,MC 

OUTPUT 

S 

PHIO«PHI (M)*CONVR 

OUTPUT 

9 

15 

WRITE  (6.3010)  M, PH 10 

OUTPUT 

10 

3010  FORMAT (*0PLANE*I4*  ANGLE  IS*F7.2*  DEGREES*) 

OUTPUT 

11 

WRITE  (6,3600)  K,Z,BNEW (M) , BZNEW (M) ,BPHNEW (M) , CNEW(M) , 

OUTPUT 

12 

1  CZNEW(M) .CPHNEW(M) 

OUTPUT 

13 

3600  FORMAT (/•  STATI0N*I5,4X*Z  IS*1PE1S,T.4X*B  1S*1PE1S.7,4X*BZ  IS* 

OUTPUT 

14 

20 

1  1PE15.7.4X*BPHI  IS*1PE15.7/7X*C  IS*1PE15,7.4X*CZ  IS*1PE15.T, 

OUTPUT 

IS 

1  4X,*CPHI  IS*1PE15.7) 

OUTPUT 

16 

WRITE  (6,3700) 

OUTPUT 

17 

3700  FORMAT  (/7X*R*12X*W*12X*U*12X*V*12X*P*10X*PT/PTO*BX*RMO*11X, 

OUTPUT 

IB 

1  *S*,12X.*M*) 

OUTPUT 

19 

25 

DO  90  N»1.NC 

OUTPUT 

20 

L*NC-N*1 

OUTPUT 

21 

IFIP(L.M) ,GT.0..AND,O(L.M> .GT.O.)  60  TO  00 

OUTPUT 

22 

AMACH*SX*XINOEF 

OUTPUT 

23 

GO  TO  BS 

OUTPUT 

24 

30 

BO  CALL  RGAS(P(L.M) ,D(L.M) ,SX) 

OUTPUT 

25 

AMACH*SORT (U (L.M) **2*V (L ,M)**2*W (L»M)**2)/AX 

OUTPUT 

26 

05  CONTINUE 

OUTPUT 

27 

PTRl«P (L,M)*(1... 5* (GAMMA-1. ) *AMACH**2) •* (GAMMA/ (GANMA-1. ) I/PTINF 

OUTPUT 

20 

WRITE  (6.3400)  R (L.M) ,W (L.M) ,U(L ,M) ,V (L.M) ,P (L.M) ,PTRL, 

OUTPUT 

29 

35 

1  D(L.M) .SX.AMACH 

OUTPUT 

30 

3400  FORMAT (1P9E13. 4) 

OUTPUT 

31 

90  CONTINUE 

OUTPUT 

32 

100  CONTINUE 

OUTPUT 

33 

RETURN 

OUTPUT 

34 

40 

ENO 

OUTPUT 

35 

1 

SUBROUTINE  RGAS  (PX.RX.SX) 

PGAS 

2 

COMMON/RGASS/AX ,HX , GX 

RGASS 

2 

C  SHORTENED  VERSION  OF  RGAS  TO  COMPUTE  ONLY  PERFECT  GAS  PROPERTIES 

PGAS 

4 

C  PX.PRESSURE  RX-OENSITY  SX-ENTROPT 

PGAS 

5 

5 

c  hx-enthalpy  ax«sound  speed 

PGAS 

6 

SX  •  ALOG(PX)  -  GX  •  ALOG(RX) 

PGAS 

7 

HX.PX/RX* ( 1 .. 1 ./ (GX-1 . ) ) 

PGAS 

8 

AX*S0RT(6X*PX/RX) 

PGAS 

9 

RETURN 

PGAS 

10 

10 

ENO 

PGAS 

11 

1 

JUMPST 

2 

SUBROUTINE  JUMPST ( OBP . DBZ . MB , NN , RR , E T AP , S I P , POM ) 

JUMPST 

3 

C  JUMPST  COMPUTES  JUMPS  CORRESPONDING  TO  DISCONTINUITIES  IN  BZ 

JUMPST 

4 

C  AND/OR  BPMI  FOR  PERFECT  GAS  ONLY. 

JUMPST 

S 

5 

C 

JUMPST 

6 

COMMON  NC,MC,K,IPRINT,P1NF,D1NF,PM!0,PI, RAO. 1. BZZ. GAMMA, M0T2.BMAX, 

CONST 

2 

1S1.S2.C1.C2, CONVR.PT INF, 

CONST 

3 

1  C(IOO) ,CZ(100) .CPNI(IOQ) ,R ( 25, 100) .0 (25. 100) ,P (25, 100) ,U(2S, 100) . 

CONST 

4 

1  V( 25,1 00 )*W( 25,100), PH 1(1001 

CONST 

5 

10 

2  .BNEW(IOO)  .CNEWUOO)  .0UMVI26)  .DUMP  (26)  , DUMU  (26)  ,OUMW  (26)  ,OUMD  (26) 

CONST 

6 

3  .BZNEW (100) .BPMNEW (100) .CZNEW(IOO) .CPHNEW ( 100) ,R0L0(26) 

CONST 

7 

4  , PTR (25. 1 00), RCLUST ( 100), PMIL( 102) 

CONST 

8 

COMMON/RGASS/AX , HX , GX 

RGASS 

2 

DATA  (INT*1B) 

JUMPST 

9 

15 

VW  ■  V (NN,MB) 

JUMPST 

10 

PW  ■  P(NN,M8) 

JUMPST 

11 

OW  *  0  (NN.MB) 

JUMPST 

12 

WW  ■  W (NN.MB) 

JUMPST 

13 

UW*U (NN.MB) 

JUMPST 

14 

2t 

CALL  RGAS(PW.OW,SWI 

JUMPST 

IS 

ASOW*AX*AX 

JUMPST 

16 

PMID-PH1 (MB) /RAO 

JUMPST 

17 

IF(IPR1NT,E0,3) WRITE  (6,3100)  NN,  PNID,  K,  Z 

JUMPST 

18 

3100  FORMAT  ( 1H0,*JUMP  IS  CALLED  FOR  AT  RADIAL  P0INT*,14,5X, 

JUMPST 

19 

25 

1  *PMI  IS*,F7.2,5X,*K  IS*. I4,5X,*Z  IS*,1PE15.6) 

JUMPST 

20 

IF ( IPRINT.EO,3) WRITE  (6,3110) 

JUMPST 

22 

3110  FORMAT (1M  ,30X,*TME  INPUT  VARIABLES  ARE  AS  FOLLOWS*) 

JUMPST 

22 

IF ( IPRINT »FQ. 31  WRITE  (6,3120)  PW,DW,UW,VW,WW,SW,ASOW 

JUMPST 

23 

A- 6 
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60 


45 


TO 


75 


SO 


es 


90 


95 


100 


105 


110 


115 


120 


125 


130 


135 


140 


5  TH£TAR*ACOS(XNMP/(XNP*XNM)) 

C0NST*M0T2-QT**2 

0SQ*QSN**2 

TMETAR*ACOS(ABS)XNMP/(XNP*XNH))) 

TMETAO»THETAR/RAO 
AMACH*SQRT (QSQ/ASQW) 

IF  (AMMSQ  .6E.  1.0)  60  TO  10 
C  SUBSONIC  CORNER  FLOW  (NO  JUMPS  IN  P.D.S.QSQ) 

QSP.QSM 
60  TO  110 

10  IF  (GSM*POM  ,LT.  0.)  60  TO  20 
C  •••  SUPERSONIC  EXPANSION  CORNER  ••• 

IF (I PRINT .£0.3) WRITE  (4.3140) 

3140  FORMAT) 1H  ,20X, ‘SUPERSONIC  EXPANSION  CORNER  WHERE*) 

!F(I PHI  NT .£0*3) WRITE  (4.3150)  TMETAO.AMACM.QT.6SM 
3150  FORMAT ( 1H  .22MTMETA.AMACM.QT.QSM  .1P4E15.5) 

CALL  R6AS(PW.0W. DUMMY) 

ANG*AS1N ( 1 ./AMACM) 

VPOWR*SQR 7 ( ASQW/ ( CONS  7-2 . *NX- ASQW  > ) 

DELoTMET AR/FLOAT ( 1NT ) 

DO  15  1*1 . INT 
CC*0. 

PW0*PW 

OWO-OW 

17  DPOAL— DW*QSQ*VPOVR 
000 AL*DPOAL/ ASQW 
CD*CC*1. 

PW> ( PWO*PW*CC*OEL#OPOAL ) /CO 
DW* (OWO.DW*CC»OEL*DOOAL ) /CO 
CALL  R6AS(PW.0W. DUMMY) 

QSQ*CONST-?.*MX 

asqw*ax*ax 

VPO VR« 1 ./SORT ( QSQ/ASQW- 1.0) 

CC*CO 

IF  (CC  ,LT.  1.5)  60  TO  17 
15  CONTINUE 

qsp*sqrt(qsq) 

60  TO  100 

C  •••  SUPERSONIC  COMPRESSION  CORNER 
20  C0STM2* (XNMP/ ( XNP*XNM) ) **2 
IF(IPRINT.E0.3IWRI7E  (4.314 Ot 

3140  FORMAT (1H  «20X**SUPERS0NIC  COMPRESSION  CORNER  WHERE*) 
IF(IPRINT.EQ.3)WRITE  (6.3150)  TMETAO. AMACM. QT.OSM 
C  •*•  (PERFECT  6AS  OBLIQUE  SMOCK  RELATIONS)  ••* 

SlNTM2*l .-C0STM2 

AM4*AMMSQ**2 

AM2*AMMSQ 

Cl  — )  (AM2*2.)/AM2»SAMMA*SINTM2) 

C3*-C0STM2/AM4 

C2* (2.«AM2*1 . ) /AM4. ( .25* (GAMMA*l. ) **2. (OAMMA-l. I /AM2) *SINTM2 
0UMM-C1/3. 

A*-C2»0UMM*C1 

SB*C3-(C2-2.*Cl**2/9.) *0UMM 
DOUMsSQRT (A/3. ) 

0DUM1*2.*00UM 

TEST  —  .5*SB/  (00UM**3) 

IF  (TEST  .GE.  -1.0)  60  TO  25 
IF(IPRINT.FQ.3)WRITE  (6.3165) 

3145  FORMAT (1M  , 20 X. ‘NORMAL  SHOCK  MODE  IS  USED*) 

PWO-PW 
T»AM2 
GO  TO  45 

25  XX*ACOS(TEST)/3. 

Xl«COS(XX) 

X2»COS)XX»2.*PI/3. ) 

X3*C0S(XX*4.*PI/3.) 

IF  (XI  .LT.  X2)  80  TO  30 

X0UM*X1 

X1*X2 

X2*X0UM 

30  IF  (XI  ,L£.  X3)  80  TO  35 
SX*X1 
60  TO  40 

35  IF  (X3  .LE.  X2)  X2«X3 
SX*X2 

40  S1NTM2*00UM|*SX-0UMM 
ANO*ASIN(SORT(SINTM2) ) 

T*AM2*SINTM2 

IF ( IPRINT.EQ.3) WRITE  (6  .  4007)  SINTH2 
4007  FORMAT  (10X.  •  SIMTH2  •  *.  1PE25.U) 

45  6A2*2.*6AMMA/ (8ANMA-1 • ) 

BO* ( 6AMMA- 1 • ) /2  » 

6E*B0«1. 


JUMPST 

53 

JUMPST 

54 

JUMPST 

55 

JUMPST 

56 

JUMPST 

57 

JUMPST 

SB 

JUMPST 

59 

JUMPST 

60 

JUMPST 

61 

JUMPST 

62 

JUMPST 

63 

JUMPST 

64 

JUMPST 

65 

JUMPST 

66 

JUMPST 

67 

JUMPST 

68 

JUMPST 

69 

JUMPST 

70 

JUMPST 

71 

JUMPST 

72 

JUMPST 

73 

JUMPST 

74 

JUMPST 

75 

JUMPST 

76 

JUMPST 

77 

JUMPST 

78 

JUMPST 

79 

JUMPST 

BO 

JUMPST 

81 

JUMPST 

82 

JUMPST 

83 

JUMPST 

84 

JUMPST 

85 

JUMPST 

86 

JUMPST 

87 

JUMPST 

88 

JUMPST 

89 

JUMPST 

90 

JUMPST 

91 

JUMPST 

92 

JUMPST 

93 

JUMPST 

94 

JUMPST 

95 

JUMPST 

96 

JUMPST 

97 

JUMPST 

98 

JUMPST 

99 

JUMPST 

100 

JUMPST 

101 

JUMPST 

102 

JUMPST 

103 

JUMPST 

104 

JUMPST 

105 

JUMPST 

106 

JUMPST 

107 

JUMPST 

108 

JUMPST 

109 

JUMPST 

110 

JUMPST 

111 

JUMPST 

112 

JUMPST 

113 

JUMPST 

114 

JUMPST 

115 

JUMPST 

116 

JUMPST 

117 

JUMPST 

118 

JUMPST 

119 

JUMPST 

120 

JUMPST 

121 

JUMPST 

122 

JUMPST 

123 

JUMPST 

124 

JUMPST 

125 

JUMPST 

126 

JUMPST 

127 

JUMPST 

128 

JUMPST 

129 

JUMPST 

130 

JUMPST 

131 

JUMPST 

132 

JUMPST 

133 

JUMPST 

134 

JUMPST 

135 

A- 7 
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G9»6C/6D  JUMPS T 

PB«Ptf*( (8A2*T>1.)/G9)  JUHPST 

DM>D»*(T/(T/89*1./GE>  1  JUMPST 

QSP«OSM«SOR  TI1.-(T-1.0)»<  8 AMMA»T  ♦ 1 . > / I T» AM2*6E*»2 ) >  JUMPST 

CALL  RGAS (PM.DM.SW)  JUMPST 

A$OM«AX*AX  JUMPST 

IE  (TEST  ,8E.  -1.0)  00  TO  100  JUMPST 

COMST-ASOM/ (6AMMA-1 .0 ) ♦ ,5»OSP»»2  JUMPST 

PK«PM«SINTM2  JUMPST 

IE  (PM  .OT.  PMO)  80  TO  TS  JUMPST 

PMbPMO  JUMPST 

TS  CALL  RGAS (PH. OB. SB)  JUMPST 

ASQH>AX*AX  JUMPST 

0SPb-S0RT(2.*(C0NST-HX>)  JUMPST 

108  D(NN.MB)«0H  JUMPST 

PINN.MB)bPM  JUMPST 

110  AAlaQT/XT  JUMPST 

AA2«QSP«POM/(DUM»XNP)  JUMPST 

V  INN. MB) >AA2* ( XNP2*SIM-XNMP*S IP) -AA1»DBZ  JUMPST 

U(NN.MB>«ABS(AA1*XT2*AA2«(XNMP-XNP2))  JUMPST 

«INN.MB)>ARS(AA1*0BP*AA2«(XNP2*ETAM-XNMP*ETAP>)  JUMPST 

IF(IPRINT.EQ.3)MRITE  (6.3170)  JUMPST 

3170  FORMATC1M  *30X.*THE  OUTPUT  VARIABLES  APE  AS  FOLLOHS*)  JUMPST 

IF(IPRINT.E0.3)«RITE(<>.3120)PM.D«.U(NN,MB>  .V  INN.  MB)  .B(NN.MB)  .SB.  JUMPST 
1  »SO»  JUMPST 

POMiANG  JUMPST 

RETURN  JUMPST 

END  JUMPST 
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APPENDIX  B.  LISTING  OF  SWINT  CHANGES 

The  modifications  extending  SWINT  to  handle  inlets  are  given  in  CDC  update 
form  in  this  section.  The  listing  contains  directive  cards  *DELETE,  *INSERT 
and  *BEFORE  which  describe  where  changes  are  to  be  inserted.  The  directive  cards 
identify  cards  in  SWINT  using  the  right  hand  column  designators  provided  in 
Appendix  B  of  Reference  2.  The  directive  cards  are  interpreted  as  follows: 

1.  *DELETE  DECK.n,DECK2.ra  .  The  cards  in  SWINT  located  between  and  including 
DECK.n  and  DECK2.m  are  deleted.  The  cards  in  the  update  listing  occurring  between 
this  *DELETE  card  and  the  next  directive  card  are  inserted  in  place  of  the 
deleted  cards. 

2.  *INSERT  DECK.m  .  The  cards  in  the  update  listing  lying  between  this 
*INSERT  card  and  the  next  directive  are  inserted  following  the  card  in  the  SWINT 
deck  with  identifier  DECK.m  . 

3.  *BEFORE  DECK.n  .  The  cards  in  the  update  listing  lying  between  this 
*BEFORE  card  and  the  next  directive  are  inserted  before  the  card  in  SWINT  with 
the  identifier  DECK.n  . 

The  listing  provided  in  this  section  also  contains  *CALL  NAME  and  *DECK  cards. 

At  the  location  in  SWINT  where  *CALL  NAME  is  added,  the  labeled  common  NAME 
should  appear.  *DECK  cards  can  be  disregarded. 

The  update  listing  contains  the  three  new  subroutines  COWL,COWLP  and  WALL2. 
The  functions  of  these  subroutines  are  analogous  to  those  of  BODY,BODYP  and  WALL, 
but  apply  to  the  cowl  surface  rather  than  the  inner  body  surface. 


uu 
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•1DCNT  COWL ADD 
•I  CC0NST.5 
•CONDECK  CCOWL 

COMMON  /CCOWL/ 

1  C(25>.  CO<2S).  CPMI (25) •  CPHI0I25).  CPHPHI(25>.  CPHPH0(2S>, 

2  CZ  (25)  *  CZO  (25)  «  CZPHK2S).  CZPHI0(25),  CZZ(25).  CZZO(25), 

3  C2MI25) •  CBM (2b) ,  C4M(25>.  C5MI25),  C7M(2S>,  ICOWL. 

A  ICOWOPTt  I JMPKTC (25) ,  IJUMPC(2S),  IJUMP1C(25>, 

5  PZCORCI25) »  SC(2S>, 

6  C3 (3) »  C4(3),  C5(3).  C7 ( 3 ) * 

7  PCY (3) t  SCY (3) *  VOCY (3) .  V?C(3> 

C 

•DELETE  CSWINT.A 

2  ASQ (20*25) •  OCT (20*25) « 

•INSERT  COOPT. 5 

*  .  ISWMODC.  ISwSMOCt  MOO 1C 

•0  SWINT.3 

1  TAPE3«51Z.  TAPE9,  TAPE  16*  TAPE1 7*512.  TAPE20,  TAPE22*512 

2  .  TAPE23*  TAPE24) 

•INSERT  SWINT.10 

•CALL  CCOWL 

•Insert  swint.3i 

IF  (ICOWL  .EO.  1)  CALL  COWL  (-11 
•INSERT  SWINT.82 

IF  (ICOWL  .EO.  1)  CALL  WALL2  (0.  M.  MP.  MM.  KN. 

S  CUP(l.NC.M).  CUP (2.NC.M)  .  CUPO.NC.MI) 

•INSERT  SWINT.92 

IF  (ICOWL  .EO.  1)  60  TO  10 

•insert  swint.ios 

IF  (ICOWL  .EO.  1)  CALL  COWL  (0) 

•INSERT  SWINT.110 

IF  (ICOWL  .EO.  1)  GO  TO  60 
•INSERT  SWINT . 163 

IF  (ICOWL  .EQ.  0)  60  TO  54 

C . ADVANCE  COWL  POINT. 

CALL  WALL2  (1.  M.  MP,  MM.  KN.  PZS,  SZS,  V2S) 

CU(l.NC.M)  a  0.5  •  (CU(l.NC.M)  «  CUP(l.NC.M)  ♦  OZ  •  PZS) 

CU(2.NC,M)  «  0.5  •  (CU (2.NC.M)  ♦  CUP (2.NC.M)  ♦  OZ  •  SZS) 

CUO.NC.N)  a  0.5  •  (CUO.NC.N)  ♦  CUP(3,NC,M)  .  OZ  *  V2S) 

54  CONTINUE 
•INSERT  SWINT. 171 

IF  (ICOWL  .CO.  1)  GO  TO  40 
•OELETE  SWINT. 193 

CALL  JUMP  (M,  1) 

•INSERT  SWINT. 196 

00  48  M  v  1  ,  MP2 

IF  (IJUMPC(M)  .EO.  0)  GO  TO  48 
IF  (M  .EO.  MP20)  GO  TO  48 
CALL  COWLP  (M.  1) 

C...0OWL  JUMP. 

CALL  JUMP  (M.  NC) 

46  CONTINUE 
•I  CORNER. 99 
•OECK  COWL 

SUBROUTINE  COWL  (JC> 

Caaaaaaaaaaaaaaaaa.aa aaaaaaaaaaaaaaaaaaaaaaaaaaaaaa .aaaa.aaaaaaaaaaaaaaa 

C 

C  THIS  ROUTINE  COMPUTES  THE  COWL  RADIUS  AND  DERIVATIVES  AT  ALL  PHI 

C  PLANES  FOR  A  GIVEN  Z.  COWLP  IS  THEN  CALLED  TO  DETERMINE 

C  NECESSARY  COWL  SHAPE  CONSTANTS  ANO  TO  CHECK  FOR  A  JUMP.  AT  THE 
C  COMPLETION  OF  COWL.  THE  COWL  RADIUS  ANO  DERIVATIVES  ARE  STORED  IN 
c  CIM) «c/ (M) .ETC.  OLD  VALUES  AT  Z  -  OZ  ARE  STOREO  IN  COIN)  CZO(M), 

C  ETC. FOR  ANY  PLANE  ON  WHICH  A  JUMP  OCCURS. NEW  AND  OLD  VALUES  ARE 

C  STORED  IN  REVERSE  OROER. 

C  NOTE  THAT  Z  IS  ASSUMED  TO  6E  INCREASING. 

C  SUBROUTINE  MUST  BE  VALID  ON  FRINGE  PLANES. 


COMMON  CCOWL 
(new) 


SUBROUTINE  SWINT 


>.  0)  GO  TO  48 


•CALL  C80DY 
•CALL  CCONST 
•CALL  CCOWL 

DO  SO  MwI.mPZ 
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..TRANSFER  NEW  VALUES  TO  OLD 
COtM)aC(M) 

CZO(M) *CZ (M) 

CPHIO(M)«CPHI (M) 

CPHPMO (M) wCPHPHI (M) 

CZPHIO (M) aCZPHI (M) 
CZZO(M)wCZZIM) 

CIM)  a  0.0 
CZ(M)  a  0.0 
CPMI (M)  a  0.0 


SUBROUTINE  COWL 
(new) 
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CPHPMI (Ml  •  0.0 
CZPHl(M)  •  0.0 
CZZ  (Ml  ■  0.0 

C***«*InSERT  DEFINITIONS  OF  C.  CZ.  CPHl.  CPHPNI.  CZPH1,  AND  CZZ  BELOB. 

C 

C 

CALL  COB LP(M.JC) 

SO  CONTINUE 
RETURN 
END 

•DECK  COhLP 

SUBROUTINE  COBlP  (N,  JC) 

C 

C  CObLP  COMPUTES  COBL  PARAMETERS  ANO  CHECKS  FOR  COBL  JUMP. 

C  INPUT  -  M  -  M-PLANE  NUMBER 

C  JC  ■  0  CHECK  FOR  JUMP  ANO  COMPUTE  COBL  PARAMETERS 

C  JC  «  -I  COMPUTES  COBL  PARAMETERS  ONLY 

C  JC  *  1  REVERSE  N£*  ANO  OLD  COBL  DESCRIPTION  AND 

C  COMPUTES  COBL  PARAMETERS 

C 

•CALL  CBOOY 
•CALL  CCONST 
•CALL  CCOML 
•CALL  CSTEPS 
•CALL  CXXYYZZ 

C  INITIALIZE  JUMP  AND  INLET  KEY 

c 

IF  (JC)  1*0.  10.  BO 
C  CHECK  FOR  A  JUMP. 

10  CONTINUE 

CZZC  ■  CZZ (M)  -  YZM(M)  •  CZPHI(M)  /  YPHIM(M) 

CZZOC  a  CZZO(M)  -  YZM(M)  •  CZPMI(M)  /  YPMIM(M) 

TEST1  ■  AMAX1  (ABS(CZZC) *  ABS(CZZOC)) 

TEST  1  MBS  (CZ  (M)  -CZO(M) I  -DZ» TEST  1 
TEST2*AMAX1 (ABS (CZPMI (Ml  I .ABS(CZPHIO(M) ) ) 

TEST2»ABS (CPMI (MI-CPHIO(M) l-DZ*TEST2 
IF(TESTl.GT.l.E-o)GO  TO  30 
IF (TEST2  .GE.  l.E-6)  GO  TO  30 
IJONPC(NI«0 
GO  TO  140 

. . JUMP  OCCURS-CMECK  FOR  OVERIOE 

30  CONTINUE 

IJUMPC (M) *1 

IF (PHI (Ml .GE.PH12J) GO  TO  HO 
if(Pmkmi.le.phiij)go  to  ao 
C.....N0  JUMP 
ICFL  •  1 
I JUMP1C (Ml «2 
IJUMPC  (MOO 
GO  TO  140 
SO  CONTINUE 
XK>CZ(MI 
CZ(MI»CZO(MI 
CZO(MI>XK 
XMCPHKNI 
CPHI(N) «CPH10(MI 
CPHI0(M1«XK 
XK-CZPMl (Ml 
CZPHI (MI«CZPMI0(MI 
CZPMIO (Ml *XK 
XKaCZZINI 
CZZIMIaCZZO(N) 

CZZOIMOXK 

XK-CPhPHl(M) 

CPHPMI (M)«CPHPM0(N) 

CPHPHO (Ml *XK 

C.... .COMPUTE  COBL  SHAPE  PARAMETERS 
140  CONTINUE 
CM*C (Ml 
CZM*CZ (Ml 

PHIZ»-YZMIMI/YPMIM(M) 

CPOB  ■  CPHIIMI  /  CM 

CPOB2  ■  CPOB  ••  2 

OUM  •  1.0  ♦  CZH  ••  2 

C22*DUM*CP0B2 

OUMl« (CPHPMI (MI/CM-CP0B2I 

DuM2* (CZZ (Ml .CZPHI (MI»PMIZI /OUM 

OUM3«CZPHI(M|/CM 

0UM4*0UM3-CZM*CP0B/CM 

C2M(N)«SONT(C22l 

C3M(MO0UMl/YPHlM(MI 


■COWLADD. 73 


-  SUBROUTINE  COWL 
(new) 


SUBROUTINE  COWLP 
(new) 


nnft 
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C4MIM) *0UM4.PMIZ*DUM1 

C5N(N)aCZPHI(H)/YPMIH(M> 

C7N(N)>0UM2*DUM 

RETURN 

END 

•INSERT  DECOOE.16 
•CALL  CCOWL 
•INSERT  OECOOE.ltS 

IF  1 ICOWL  .EO.  01  60  TO  2000 


SUBROUTINE  COWLP 
inew) 


...KCOOE  COWL  POINT. 


U3  ■  CVI3.NC.N) 

SCN  «  CV (2.NC.M) 

SC  (Ml  «  SCN 
T1  ■  CZ(M) 

T2  »  CPHIINI  /  C (Ml 
T3  •  1.0  «  T2  *•  2 
PM  ■  EXP  (CV(l.NC.M)) 

IF  (CV < 1 «NC*M)  »6T •  -600.0  .AND.  CV(l.NC.M)  .LI.  700.01  GO  TO  2002 
WRITE  (6  t  34571  M.  PM 
CALL  SAVE 

2002  CONTINUE 
P(NC«MI  -  PH 

CALL  R6AS  (PM.  DM.  SCM.  6) 

D(NC.M)  ■  OH 

ASO(NC.M)  •  AX  •  AX 

(ISO  *  MOT 2  -  2.0  •  HX 

CDUMP  >  osa  •  T3  -  U3  •  U3 

IF  (COUMP  .OE.  0.01  GO  TO  2003 

CALL  DHPSQRT  (6HDEC00E.  5.  Z.  K.  M.  NC.  COUMP) 

2003  CONTINUE 

WM  «  SORT  (COUMP)  /  C2-(M) 

W(NC.M)  -  WM 

V (NC t M)  •  (U3  -  T2  •  TI  •  AM)  T3 
U(NC.M)  «  Tl  «  WM  .  T2  •  V(NC.M) 

GO  TO  9 
2006  CONTINUE 
•INSERT  DECODE .233 

IF  (ICOWL  .EO.  1)  GO  TO  100 
•INSERT  0EC00E.236 
100  CONTINUE 
•INSERT  DECODE. 246 

3457  FORMAT  (1H1.  •  IN  SUBROUTINE  DECODE  THE  LOG  OF  PRESSURE  ON  PLANE*. 

S  14,  •  ON  THE  COWL  IS  •.  1PE1S.6,  5X.  - -  STOP  — •) 

•INSERT  EVAL.1S 
•CALL  CCOWL 
•INSERT  EVAL.29 

IF  (ICOWL  .EO.  0)  GO  TO  15 
C3(KF>  ■  C3M (M) 

C4(KF)  ■  C4M(M) 

C5(KF>  •  C5M(MI 
C7(KF)  ■  C7M(N) 

IS  CONTINUE 
•INSERT  EVAL.101 

IF  (N  .LT.  NC  .OR.  ICOWL  .EO.  0)  GO  TO  26 
PCY(KF)  «  PNN 
VOCY(KF)  *  VNN  /  WNM 
SCY(KF)  ■  SC (N) 

V2C(KF>  ■  VNN  ♦  CPHI(H)  •  UNM  /  C(M> 

26  CONTINUE 
•INSERT  FIELD. 15 
•CALL  CCOWL 
•CALL  COECOOE 
•INSERT  FIELD. 74 
6A*6A2* .5 

PTM  •  (1.0  *  GD  •  ACH**2.0)  ••  (-6A) 

IF ( ICOWL.EO.O) GO  TO  40 
PTRTOT  »  0.0 
AREAT  *  0.0 

C... ..PRESSURE  RECOVERY  CALCULATED 
C.....PTR  IS  LOCAL  TOTAL  PRESSURE  RATIO 

C..... PRESSURE  RECOVERY  ■  SUM  (PTR  •  LOCAL  AFFECTED  AREA)  t  TOTAL  AREA 
00  111  N  •  2.  MPl 
INI  ■  1 
INC  •  1 

IF  ( (1-IDYAW)*N  .EO.  2)  INI  ■  0 
IF  ( (1-IDYAW)«N  .EO.  MPl)  INC  ■  0 
00  110  N  ■  1.  NC 

AX  ■  SORT (GAMMA  •  P (N.M) /D (N,M>> 

AMACH  ■  SORT (U(N,M>  **2.0  •  V(N,M)*«2.0  *  W(N,M)«*2.0>  /  AX 
PTR  *  PIN. Ml  •  PTN  •  (1.0  ♦  GO  •  ANACM**2I ••OA 
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I”  SUBROUTINE  EVAL 


f-  SUBROUTINE  FIELD 


B-4 


NSWC  TR  83-428 


; 


INI  a  1 
INC  a  1 

IF  (N  .CO.  1)  INI  ■  0 
IF  (N  .CO.  NC>  INC  •  0 

AREA  ■  0.25  •  (PHI(N*IMC)-PHI(N-IN1))  •  I ( O.S* (R (N.m) * 

I  H(N.INC.N))) **2.0  -  (0.5*(R(N-IN1,M)«R(N.M) ) >**2.0> 

AREAT  ■  AREAT  *  AREA 
PTRTOT  a  PTRTOT  ♦  PTR  •  AREA 

110  CONTINUE 

111  CONTINUE 

PR  a  PTRTOT  /  AREAT 
MRITE  16.3005)  PR 
40  CONTINUE 
•INSERT  FIELD. 84 
PTRaXINOEF 
•I  FIELD. 88 

PTR-P  CL.N) *PTN* ( 1 . *60«AMACH**2> ••6A 
•OELETE  FIELD. 99 
•OELCTE  FIEL0.100 

MRITE (6, 3400 IL.R(L.M) .M(L.M) .UIL.M) , V <L .M) .P ( L .M) , PTR , 

I  0(L.H),SX.AMACH,TN1,TZ1,ISX 
•INSERT  FIELD.114 
PTRaXINOEF 
•INSERT  FIEL0.118 

PTR«P  <L .MM) »PTM« 1 1 . ♦G0*ANACH««2> ••6A 
•DELETE  F1EL0.120 
•DELETE  FIELD. 121 

MR1T£(6.3400)L,R(L,NM) .M(L.MM) .U(L.MM) .V(L.MM) .P(L.MM) • 

)  PTR.O (L.MM) .SX* AMACH.THR (lCF.D.THZdCF.Ll.ISOCF.L) 

•INSERT  FIELD. 129 

3005  FORMAT  (*0*.*PRESSURE  RECOVERY  PTAVR/PT1NF  •  •.F10.5) 

•OELETE  FIELD. 131 

3400  FORMAT (*  •.12.1X,3(1P£L1«4,1X) .1P8E11.4.1X.14) 

•DELETE  FIELD. 136 
•OELETE  FIELD. 13T 

3700  FORMAT  1*0  N  •.6X.1HR.11X.1HM.11X.1HU.11X.1HV.10X.1HP.TX.BHPT/PTINF. 
1  6X.3HRHO.9X.1HS.10X.1HM.10X.2HTN.9X.2HTZ.5X.2HIS) 

•INSERT  FINA0.1S 
•CALL  CCOML 
•INSERT  FINAD.8T 

IF  (ICOML  .EQ.  1)  CALL  MALL2  *0.  M2,  HP.  MM,  KN« 

S  CUPd.NC.M2),  CUP (2.NC.M2)  «  CUP (3.NC.M2) ) 

•INSERT  FINAO.llT 

IF  (ICOML  .CO.  I)  80  TO  35 
•INSERT  FINA0.199 

IF  (ICOML  .EO.  0)  60  TO  61 

CALL  MALL2  (1.  M2.  MP.  MM,  KN.  PS,  SS.  VS) 

CUU.NC.M2)  a  0.5  •  (Cud.NC.M2)  .  CUP(1.NC,M2>  ♦  PS  *  DZ) 

CU(2,NC.M2I  a  o.S  •  (CU(2,NC.M2)  *  CUP(2,NC.M2)  »  SS  •  DZ) 

CUO.NC.M2)  a  0.5  •  (CUO.NC.M2)  .  CUPd.NC.M2)  *  VS  •  OZ) 

61  CONTINUE 
•INSERT  F1NA0.231 

IF  (ICOML  .EO.  1)  60  TO  TO 
•INSERT  FRINGE. 10 
•CALL  CCOML 
•BEFORE  FRIN6E.19 

IF  (ICOML  .EO.  0)  60  TO  93 
CPHPHI  (Ml )  •  CPHPHKM2) 

CZZ(Ml)  a  CZZ (M2) 

CZPMI  (Ml )  «  CZPMKM2) 

SC(M1)  a  SC (M2) 

93  CONTINUE 
•INSERT  INIT.6 
•CALL  CCOML 
•INSERT  IN1T.83 

IF  (ICOML  .EO.  0)  60  TO  42 
SC(M)  a  SFF 

CU(l.NC.M)  a  AL06  (P(NC.M) ) 

CU12.NC.M)  a  SOM) 

CUO.NC.M)  a  V  (NC.M)  .  U(NC.M)  •  CPHI  (M)  /  C(M) 

60  TO  43 

42  CONTINUE 
•INSERT  INIT.86 

43  CONTINUE 
•INSERT  1NIT.99 

IJUMPC(M)  a  0 

IF  IICOMOPT  .EO.  1)  60  TO  33 
IJUMPIC(M)  a  0 
IJMPKTC(M)  a  0 
00  TO  34 
33  CONTINUE 

lUUMPlC(M)  •  4 
I JMPKTC (M)  a  1 
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34  CONTINUE 
•INSERT  INLET. 14 
•CALL  CCONL 
•DELETE  INLET. 70 

WRITE  (6.  1000) 

1000  FORMAT  '•  INLET  SUBROUTINE  CALLED  •) 

CALL  JUMP  (M.  1) 

•INSERT  INTE6.14 
•CALL  CCOWL 
•INSERT  INTE6.112 
C 

IF  (1COWL  ,E0.  0)  60  TO  1600 
...INTEGRATE  THE  COWL  PRESSURES. 

Kl-1 

DO  1325  1-1.6 
SSFlD-0. 

SUMJU.  11-0. 

SUMJI 1,2) >0. 

132S  CONTINUE 

IF  IIOYAW  .EQ.  II  60  TO  1300 
00*  SIMPSON-S  RULE  FOR  SYMMETRY  CASE  (PMIO-180)  ••• 

DO  1200  M-2.MC 
CM  •  C (M) 

CPH18  ■  CPHI(M)  /  CM 

DUM  c  >2.0  •  (P(NC.M)  -  PINF)  •  CM  •  PHlO  •  TG4M<M> 

SINP  a  SINPHIIMI 
COSP  -  COSPHIIM) 

SF  1 1  > -DUM<  (C0SP*CPHIB*S1NP) 

SF (3) -OUM*CZ (Ml 
DUM1-SF (3)*CM 
SF (2) *DUM1*C0SP 
IF (M. ME. 3)60  TO  112S 
SF 1  €  1 )  a  SUM J (1.1) 

SF1I2)  a  SUMJI 2,1) 

SF1 (3)  a  SUMJI3.il 
1126  DO  1150  1-1,3 
1150  SUMJII, K1)-SUMJ(I,K1).SF(I> 

K1-3-K1 
1200  CONTINUE 

SF2II)-2.»IPINC.MPI|-PINFI»C(MP1)«PHI0»TG4M(MP1) 

SF2I3) --SF2I 1)*CZ (MP1 I 
SF2 (21 a-SF2 (3) *C(MPl I 
DO  1250  1-1,3 

FID -0Y03* (4.-SUMJ 1 1 ,2) «2.*SUMJ (1.11) -SSF  < I) *F  < I) 

IF  (K1  .10.  1)  GO  TO  1225 

.....EVEN  NUMBER  OF  POINTS-USE  SIMPSONS  RULE  FROM  M-2  TO  MC. 

APPROXIMATE  THE  LAST  INTERVAL  USING  TRAPEZOIDAL  INTEGRATION. 
F(II-F(I).DY03»(.5«SF(II-SF1UI»1.5»SF2<I)> 

60  TO  1250 

1226  F 1 1 1 -F 1 1 ) »OYD3» (SF2 1 1 ) -SFI ( I ) ) 

C.....OOD  NUMBER  OF  POINTS-  APPLY  SIMPSONS  RULE  DIRECTLY. 

1250  CONTINUE 

60  TO  1600 

C  •••  SIMPSON-S  RULE  FOR  NON-SYMMETRIC  CASE  (PHIO-360)  ••• 

1300  CONTINUE 
Kl-1 

DO  1500  M-2.MP 1 
CM  a  C1MI 

CPMIB  -  CPHIIM)  /  CM 

OUM— (  P  ( NC  .  M ) -P I NF  )  «CM»PM  1 0*  T64M I M ) 

SINP  a  SINPHIIMI 
COSP  a  COSPHIIM) 

SF I 1 I aDUM* I COSP .CPH IB* S 1 NP ) 

SF (31  a  OUM  •  CZ (Ml 

OUMI  a  SFI3I  •  CM 

SF (2)  a  OUMI  •  COSP 

SFI*)  a  OUMI  •  SINP 

SFI5) -OUM* (CPHIB*COSP-SINP) 

SF  1 6 1  aDUM*CM*CPM  I B 
IFIM.NE.3I60  TO  13T5 
DO  1350  1-1,6 
SF1II)  a  SUMJII, 1) 

SF2II)  a  SF 1 1 ) 

1350  CONTINUE 

1376  DO  1400  I>1,6 

1400  SUMJII, Kl)aSUMJ(I,Kl)*SF(I) 

K1 -3-111 
1500  CONTINUE 

DO  1525  1-1,6 

Fill «0Y03* I4.*SUMJ( 1,2) *2,*SUMJ( 1,1) )»SSF 1 1 1 ,F 1 1 ) 
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C.....000  NUMBER  OF  POINTS-APPLY  SIMPSONS  RULE  DIRECTLY. 

IF  (Ml  .EQ.  1)  60  TO  1525 

. . EVEN  NUMBER  OF  POINTS-APPLY  SIMPSON  RULE  FROM  M»2  TO  MP1  AND 

C  EXTEND  TO  Ma3.  USE  TRAPEZOIDAL  RULE  TO  SUBTRACT  Ma2  TO  3  INTERVAL. 
F i I I «F( IT ♦0Y03* (A.*SF1 ( I ) -SF2 ( I) -2.*SF ( I) ) 

1526  CONTINUE 
C 

1606  CONTINUE 
•DELETE  JUMP. 2 

SUBROUTINE  JUMP  (MB.  NT 
•INSERT  JUMP. 16 
•CALL  CCOML 

•DELETE  JUMP .2 A. JUMP. 26 
SO  >  SV(MB) 

POM  a  1.0 
ETAP  a  BZ (MB) 

SIP  a  BPHI (MB)  /  B(MB) 

IF  (N  .EO.  II  60  TO  AO 
C...OOWL  SURFACE. 

OBP  a  CPNIO(Mfi)  /  CO  (MB I  -  CPHKMST  /  C(MB) 

OBZ  a  CZO(MB)  -  CZ(MB) 

SO  a  SC (MB) 

IF  (N  .EO.  I)  WRITE ( JJJJJ.3130 1  6 (MB) . BZ (MB) .BPHI (MB) tOBP.OBZ.HOT2 
IF  (N  .6T.  1)  WRITE ( JJJJJ.3131 )C (MB) .CZ (MB) .CPH1 (MB) .DBP.DBZ.H0T2 
POM  a  .1.0 
ETAP  a  CZ  (MB) 

SIP  a  CPHI(MB)  /  C (MB) 

AS  CONTINUE 

UW  a  U(N.MB) 

VW  a  V(N.MB) 

WN  a  W(N.MB) 

PM  a  P(N.MB) 

DW  a  DIN. MB) 

ASQW  a  A SO  IN. MB) 

•OELETE  JUMP. 32 
•DELETE  JUMP.A9 

IF  IN  .EO.  1)  I JUMP 1 (MB)  a  2 
IF  IN  «6T •  1)  IJUMPIC(MB)  a  2 
•OELETE  JUMP. 55 

IF  IN  .EQ.  1 >  I JUNPJ I MB)  a  2 
IF  (N  .GT.  II  I JUMP1C (MB)  a  2 
•OELETE  JUMP.5T 
10  CONTINUE 

IF  (OSM*POM  .LT.  0.01  GO  TO  20 
•OELETE  JUMP. 62 
ICFL  ■  1 

IF  IN  .6T.  1)  60  TO  15 
IJUHP1 (MB)  a  3 
I JMPKT (MB)  a  1 
60  TO  100 
16  CONTINUE 

IJUNPIC(MB)  a  3 
IJMPKTC(MB)  a  i 
•OELETE  JUMP. 66 .JUMP. 68 

IF  (N  .GT .  1)  GO  TO  22 
IJUMPl(MB)  a  A 
I JMPKT (MB)  a  1 
60  TO  23 

22  CONTINUE 
IJUMPIC(MB)  ■  A 
I JMPKTC (MB)  ■  1 

23  CONTINUE 

CALL  COMP  (THETAR.  AMACH.  PW»  DW.  SO.  QSM,  OSP,  OT .  ASOW.  XUK) 

100  CONTINUE 

DIN, MB)  «  DW 
PIN, MB)  a  PW 
ASOIN.MB)  a  ASOW 
•OELETE  JUMP. 70. JUMP. 72 

IF  (N  .EO.  II  SB (MB)  a  SO 
IF  (N  «6T .  1)  SC ( MB )  ■  SO 

V(N.MB)  •  AA2  •  (XNP2  •  SIM  -  XNMP  •  SIP)  -  AA1  •  DBZ 

U(N.MB)  ■  AAt  •  XT 2  ♦  AA2  •  (XNMP  -  XNP2I 

W(N.MB)  a  AA1  •  DBP  ♦  AA2  •  (XNP2  •  ETAM  -  XNMP  •  ETAP) 

•OELETE  JUMP. 7A, JUMP. 77 

WRITE  (JJJJJ.31201  PW.DW.U(N.MB) .VIN.MBI .WIN. MB) ,S0, ASQ (N.MB) 

CUIl, N, MB)  a  AL06  (P(N.MBI) 

CUI2.N.MB)  •  SO 

CUO.N.MB)  ■  V (N.MB)  «  SIP  •  U(N.MB) 

•INSERT  JUMP. BA 

3131  FORMAT  (1H  .22MC.CZ.CPHI .DBP.DBZ.H0T2, 1 POE  15.5) 
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•INSERT  OUT. • 

•CALL  CCOHL 
•INSERT  OUT. 61 

HRITE  1231  Z»  NP2N1.  (P8(M),  N-1.MP2M1) 

•INSIRT  OUT.TI 

IF  (ICOHL  .EQ.  0)  60  TO  650 
C 

C...1MIS  SECTION  OUTPUTS  PRESSURES  AT  THE  COHL. 

C 

REMIND  16 
NPTS  >  0 
Ml  -  1 
NCNX  •  MC 
N1P14  ■  HI  ♦  16 
T2S  CONTINUE 

NPTS  -  NPTS  *  1 
730  CONTINUE 
NCS  *  NC 

REAO  <161  NC.MC.ATTA, YAH,ACH,GAMMA,PINF,DINF .PH10.K.Z. 

1  (OUM.  1-1.12).  (PHI (Ml.  DUM,  OUM.  DUM.  M-1.MP2M1). 

2  ( (DUH.  DUM,  D UN,  OUM.  PBIM).  DUM,  H-1.MP2M1).  N-l.NO 
IF  (EOF  (16)1  640,  750 

750  CONTINUE 

M2  *  NINO  (MC.  H1P14) 

IF  (H2  ,LT.  Ml)  60  TO  730 
IF  (NC  .NE.  MCSI  NPTS  >  1 
IF  (HOD  INPTS-1 ,  S3)  .NE.  0)  60  TO  775 
HRITE  (6,  3000)  ACH,  ATTA.  YAH.  20 
IF  (IPC1D  .EQ.  0)  UNITE  (6  ,  3520) 

IF  (IPCIO  .NE.  0)  HRITE  (6  .  35261 
DO  760  H  ■  Ml  .  M2 

PHI (M)  >  PHI (Ml  /  NAD 
760  CONTINUE 

HRITE  (6  .  3030)  (PHI (HI ,  M-M1.M2) 

HRITE  (6  .  30401 
776  CONTINUE 

00  762  M  ■  Ml  .  M2 

PB INI  ■  PBIM)  /  PINF 
762  CONTINUE 

zz  ■  z  •  zo 

IF  (IPCIO  .EQ.  1)  60  TO  800 

HRITE  16  .  3050)  ZZ.  (PBIM).  M-H1.M2) 

60  TO  825 
800  CONTINUE 

DO  810  M  ■  Ml  ,  M2 

PBIM)  ■  (PBIM)  -  1.0)  /  CONTI 
810  CONTINUE 

HRITE  (6  .  30551  ZZ.  (PBIM).  M-M1.N2) 

626  CONTINUE 

IF  (Ml  ,6T.  1)  60  TO  725 
NCNX  »  MAXO  (NC,  MCMX) 

HRITE  124)  Z,  HP2H1  •  (PBIM).  MM.MP2M1) 

60  TO  725 
040  CONTINUE 

Ml  ■  Ml  •  15 

IF  (Ml  ,6T .  MCMX)  60  TO  850 

M1P14  •  Ml  ♦  14 

NPTS  ■  0 

REMIND  16 

60  TO  725 

C4* *••***•*••• •••**••••••*••••••*•**•*•••*••**••••••*••••••••••*•••••*•• 

650  CONTINUE 
•INSERT  OUT. 262 

3520  FORMAT  (1H0,35X.«C  OHL  PRESSURE  RATIO*) 

3526  FORMAT  (1H0,30X,*C  OHL  PRESSURE  •> 

♦  *C  0  E  F  F  I  C  I  E  N  !•) 

•BEFORE  PR1NTST.10 
•CALL  CCOML 
•DELETE  PRINTST.26 

1  ICOHL.  NSFO,  NS60 
•INSERT  PR1NTST  33 

IFIICOML.EQ.DHRITE  (6  ,  5023)  ISHSMOC,  ISHMODC,  M0D1C,  ICOHOPT 
•INSERT  PRlNTST .94 

6  M.«  ICOHL  ■  •.I10.I0X, •(OUTER  80UN0ARY  DEFINITION.  1«HALL.0»* 

S  •SHOCK) •./, 

•INSERT  PRlNTST. 157 
C 

5023  FORMAT  (//✓.  5X,  4M««»»»  •  COHL  OPTIONS*.  4M»»»»,  //, 

1  OX,*  ISHSMOC  ■  •• 16. 10X.«( ISHSMOC  ■  1  -ENTROPY  EXTRAPOLATION,  ■ 
SO— ST  ANOARO I • . / , 

2  OX,*  ISUMOOC  ■  *, 16, 1 0X,« (FORM  OF  8OUN0ARY  CONDITIONS-  0  ■  14A. 
S1SA,  3  >  14C.15C)*./, 

3  OX,*  MOOIC  «  *,I8,10X** (ORDER  OF  ACCURACY—  0*1ST  ORDER,  1*2ND 
SORDER  UNTIL  COHL  DISCONTINUITY  ENCOUNTERED) •«/* 
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4  9X**  ICOMOPT  »  **  16* l  OX**  (COML  NORMAL  DERIVATIVE  CONTROL-  0>* 

5  * STANDARD* (^MODIFIED  FOR  4  STEPS)*) 

•INSERT  RADIUS.? 

•CALL  CCOML 
•DELETE  RADIUS. 20 

IF  ( I  COWL  .EQ.  0)  CC  ■  CUU.NC.M)  -  8(H) 

IF  (ICOHL  .EQ.  1)  CC  «  C(M>  •  B(M) 

•INSERT  READIN.10 
•CALL  CCOML 
•INSERT  READIN.34 

6  .  ICOHL *  ICOMOPT*  ISMMOOC*  ISMSMOC.  M001C 
•INSERT  REAOIn.R) 

ICOHL  ■  0 
ICOMOPT  ■  0 
ISMMOOC  *  3 
ISMSMOC  •  0 
M001C  «  1 

•INSERT  REA0IN.15T 
KFAC1  >  1 
•INSERT  REZONE. 16 
•CALL  CCOML 
•INSERT  REZONE. 68 

IF  (ICOHL  .EQ.  1)  CALL  COML  (-1) 

IF  (ICOHL  .EQ.  1)  60  TO  65 
•INSERT  REZ0NE.T8 
66  CONTINUE 
•INSERT  SAVE. 10 
•CALL  CCOML 
•INSERT  SHOCK. 14 
•CALL  CCOML 
•INSERT  TRANF.1T 
•CALL  CCOML 
•INSERT  STEPS. 20 

C.....AT  LAST  STEP*  MAKE  SURE  Z>ZEND. 

IF  (Z*DZ  ,LE.  ZEND)  60  TO  204 
DZ  ■  ZENO  -  Z 
Z  ■  ZENO 

IF  (OZ  .LT.  l.E-4)  OZ  «  l.E-4 
60  TO  205 
204  CONTINUE 
•I  MALL. 140 
•DECK  MALL2 

SUBROUTINE  MALL2(L*M* JR* JL.JSG.PZ.SZ*  V2Z) 

C 

C  MALL2  COMPUTES  PREDICTED  OR  CORRECTED  Z  DERIVATIVES  OF 

C  P.  V2.  AND  S( ENTROPY)  US1N6  CHARACTERISTIC  COMP.  RELS. 

C  V2  IS  VEL.  COMP.  TANGENT  TO  COML  V2«V* (CPHI/B)«U 

C  V2C ( J) •  SCY ( J) •  VOCYIJ).  AND  PCYIJ)  ARE  COML  VALUES  OF 

C  V2«  S.  V/M  ANO  P  RESPECTIVELY.  CONTROL  INTEGERS  IN  ARGUMENT  ARE 
C  L  .  *0  FOR  PREDICTOR  AND  1  FOR  CORRECTOR 

C  M  ,  Y  PLANE 

C  JR.  STORAGE  LOCATION  OF  RIGHT  SIDE  DIFFERENCE  QUANTITIES 

C  JL.  STORAGE  LOCATION  OF  LEFT  SIDE  DIFFERENCE  QUANTITIES. 

C  JS6.  STORAGE  LOCATION  OF  M  PLANE  FOR  DIFFERENCED  QUANTITIES. 

C  JR  ANO  JL  ARE  LINE  IOENT.  INDEXES  FOR  Y  OIFFS. 

C  JS6«1.2.3  LINE  1N0EX  FOR  TRANG  ANO  COML  PARAMETERS 

C  IF-1.2  LINE  INDEX  FOR  TKANF  PARAMETERS 

C  THIS  VERSION  OF  MALL2  CONTAINS  SEVERAL  OPTIONS  FOR  COML  B.C. 

C  ISMSMOC  NE  0  MEANS  COML  ENTROPY  EXTRAPOLATION 

C  M001C  *  1  MEANS  SECOND  ORDER  ACCURACY 

C  ISMMOOC  *  0  MEANS  MOO  0  FOR  COML  B.C. 

C  >3  MEANS  MOO  3  FOR  COML  b.C. 

C  THIS  ROUTINE  CONTAINS  SPECIAL  FEATURES  AFTER  A  JUMP 

C  I JUMP IC (M)  *  0  MEANS  NO  JUMP  ON  LINE 

C  I JUMP 1 C ( M )  NE  0  MEANS  JUMP  HAS  BEEN  CALLED  (SEE  JUMP) 

C  IJUMPIC(M)  ■  2  MEANS  NO  SECOND  ORDER  ACCURACY 

C  AND  NO  ENTROPY  EXTRAP.  IF  A  COMPRESSION  JUMP 

C  ANO  MOO  0  FOR  COML  B.C. 

C 

•CALL  CSmINT 
•CALL  CCONST 
•CALL  CCOML 
•CALL  COECOOE 
•CALL  COOPT 
•CALL  CEVAL 
•CALL  CSHOCK 
•CALL  CSTEPS 
•CALL  CTRANF 
•CALL  CTRANG 
•CALL  CMALL 
•CALL  CXXYYZZ 
C 

DIMENSION  0C6YU)  .ICONT(IOO) 


J-  SUBROUTINE  PRINTST 
-SUBROUTINE  RADIUS 


L SUBROUTINE  READIN 


H SUBROUTINE  RE70NE 


SUBROUTINE  SAVE 
|- SUBROUTINE  SHOCK 
h  SUBROUTINE  TRANF 


SUBROUTINE  STEPS 


U  SUBROUTINE  WALL2 
(new) 
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OAT*  ICONT/100*0/ 

KMOOIC  ■  MOO 1C 
KSWMOD  ■  ISNMOOC 
CN  ■  C(M) 

C2M  «  CZ(M) 

CPMOB  a  CPMKN)  /  CM 

YPHIJaYPHKJSGI 

YZJ-YZIJSB) 

XRVaXR(NC.JSG) 

UN  a  U(NC,M) 

NN  ■  N(NC.M) 

PN  »  PCMC.M1 
VN  >  V(NC.M) 

ON  »  O(NC.M) 

ASON  a  ASO(NC.N) 

VOR  *  VN  /  CM 

VON  a  VN  /  NN 

YPOR  a  YPMIJ  /  CM 

B8aWN*YZJ.YPMlJ*V0R 

PX  ■  (PN  -  PINA. Ml)  •  OOX 

DNN  «  ON  •  NN 

CTA  *  NN  ••  2  -  ASON 

COUMP. (ETA* ( 1 .«CPMOe**2l ♦ ( NN*CZM) /ASON 

IF  (COUMP  .6E«  0.)  GO  TO  1001 

CALL  OMPSOMT (5HMALL2. 1 .Z.K.M.NC. COUMP) 

1001  BETA  a  -SORT (COUMP) 

ALAM  a  ASON  •  (BETA  -  CZM)  /  ETA 
DUMA  a  VN  ♦  CPHOB  •  UN 
AAX  a  -UNOR(NA.JSG)  •  DOX 
DSY*(SCY ( JR) -$CY ( JL)  )  *ODT 
V2Ya  ( vac  ( JR > -V2C  ( JL > » *OOY 
IF  (IJUNPIC(M)  .EG.  01  GO  TO  20 
IF  ( IJUMPIC(M)  .tO.  2)  60  TO  15 
IF  (L  .EQ.  1)  GO  TO  10 
IF  ( 1 JUMP1C (Ml  .ME.  31  GO  TO  210 
IF  ( I JMPKTC (Ml  .GT.  NJNPKT)  GO  TO  230 
IFAC  a  NJMPKT 
60  TO  215 
210  CONTINUE 

IF  ( I JMPKTC (Ml  .GT.  NJMKTCI  GO  TO  230 
IFAC  a  NJNKTC 

IF  (ICFL*KCFL  .NE.  0)  GO  TO  215 
I JMPKTC (Ml  a  I JMPKTC (Ml  .  KFAC 
GO  TO  10 
215  CONTINUE 

I JMPKTC (Ml  a  I JMPKTC (Ml  ♦  1 
GO  TO  10 
230  CONTINUE 

IJUMP1CIM1  a  2 
I JMPKTC (Ml  a  0 
GO  TO  15 

10  FACaFLOAT ( I JMPKTC (Ml -1 1/FLOAT ( IFAC) 

PX  a  FAC  •  PX 
AAX  a  FAC  •  AAX 

HRITE  (LLLLL.2002)  Z.  M.  L,  FAC 
16  CONTINUE 
MOOIC  a  o 
ISNMOOC  >  0 
1CONT (Ml alCONT (M)  ♦  1 

IF  ( ICONT (Ml  .EG.  II  WRITE  (LLLLL.2000)  Z.  MOOIC.  ISNMOOC.  M 
20  CONTINUE 

IF  (ISNMOOC  .EG.  31  50  TO  25 

OVONY  a  (VOCYIJRI  -  VOCYIJLIl  •  DOT 

PY  •  (PCY ( JRI  -  PCY ( JL) 1  •  OOY 

OUNM  a  (CZM  /  CM  ♦  YPOR  •  OVONY)  •  NN 

DUM  a  BB  •  NN  /  ASON  -  YZJ 

PZ23»(ALAM*0UM»CZM*YZJ.YP0R*CPM0BI *PY 

1  -ONN*  (BB*  (C5  ( JSG)  *V0N*C3  ( JSG) ) 

2  —At tp* (OUMM. ( TF6 (NC» JSG) — T66 ( JSG) I *NN*VOR*  TFT (NC . JSG) 1 1 
V2Z23  (BB* (V2Y-UM*C3 ( JSG) l«YPOR*PY/DNI /NN 

GO  TO  50 
26  TG5JaT65 (JSG) 

00  AO  I ■ 1 . A 

OCGY(Z)a(CO(I.NC.JR)*OETY(NC.JR)-CG(I.NC.JL)*DETY(NC.JL))*DOY 
AO  CONTINUE 

CE1*0N* (UN/CM.TG5J*BB*TF6(NC. JSG) *WN.TF7 (NC. JSG) *VOR) 

0UM1  a  (T65J  •  YZJ  ♦  TOO l JSG) I  *  PN 
DUM2  a  (TG5J  •  YPHIJ  *  CPMOB)  •  PN  /  CM 
GSG»NN**2.VN**2.UN**2 

OUM3*UN*OCGY (31 -OSG*OCGY ( 1 1 »NH* (OCGY (2) *0UM1 ) »VN* (DCGT (A) .0UM2) 
C  Of*  TME  EQUATION  FOR  XK1  IS  VALIO  FOR  PERFECT  GAS  ONLY  ••• 
XK1*-0M/(ASGN*G8I 


SUBROUTINE  WALL 2 


B-10 
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PZ23*ALAM*WW* (C£l «2.*DCGY ( 1 I «XKl*DUM3/0* ) 

1  ♦ (CZM-ALAM) • (DC6Y (2) *0UM1 ) -OCGY  <3) ♦CPHO0* (DCGY (A) *0UH2> 

V2Z23*  (CPHOB*OCGY  (3)-0UM4*0C6Y  ( 1  >  *DC6Y  ( 4 )  *0UM2 1  /DWV 
59  PZ»ALAN*XRW*PX-(ONW*(ALAM*AAA-WW*C7(JS6l 
1  -VW*C4 ( JS6) *DUM4* VOW/CM) *PZ23)/BETA 

V2Z»UW*C4  <  JSG) - VOR*CZM-V2Z23 
IF  (MOD 1C  .N£.  1)  60  TO  90 
IF  (L  .CO.  1>  60  TO  60 

PXX  ■  (-2.0  •  P(NA.M)  •  P (NC2.M)  ♦  PA)  •  OOX 
AAXX  •  (-2.0  •  UNOP (MA. JSG)  ♦  UNOR (NC2.USG) I  •  OOX 
PZCORC(M) «ALAN* (XRW*PXX-OWW* AAXX/BETA I 
60  TO  90 

St  PZ«PZ«PZCORC(M) 

90  IF  (M  ,6T •  ISWSNOC)  60  TO  100 

CALL  R6AS(P(NC2*M) .0(NC2.M) .SW3.4) 

CALL  RGAS(P(NA«M) .D(NA.M) .SW2.4) 

SZ*2.*SM2-SW3 

IF  (SW2.LT.  Sw3)SZ«.5*(SW2*SW3) 

CUP (2.NC.M) «SZ*FLOaT (L> 

CU(2.NC.M)*SZ 

SZ-0. 

GO  TO  125 
100  SZ*-BB*DSV/WW 
126  CONTINUE 
PZ-PZ/PW 
MOD 1C  «  KMOOIC 
ISMMODC  »  KSWMOD 
IF (IJUMPC(M) .£0.0)60  TO  110 
IF (L.E0.D60  TO  110 
PZ»0. 

WRITE  (LLLLL.2001)  Z.  M 
110  CONTINUE 
RETURN 

2000  FORMAT  (SX.  •  FROM  WALL2—  AT  Z  *  (.FIS.?.*  MOOIC  AND  ISWMODC  ARE  PE 
SRMANENTLY  SET  TO*. 214. •  ON  PLANE  ••  14) 

2003  FORMAT  (SX.  •  FROM  WALL2—  IN  PREOICTOR  STEP  AT  2*  *,  F10.S. 

1  •  PZ  IS  SET  TO  0.0  ON  PLANE*.  IS) 

2002  FORMAT  (5X.  •  FROM  WALL2—  AT  2*  •»  FlS.7.  •  AND  ON  PLANE*.  15* 

1  •  L«  *.  IS.  •  X  DERIVATIVES  ARE  SCALED  BY  *.  F10.5) 

END 
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APPENDIX  C.  USER  INSTRUCTIONS  FOR  APPLICATION  OF  SWINT 
TO  AN  INLET  CONFIGURATION 

The  first  step  in  calculating  an  inlet  configuration  is  to  determine  the 
flow  field  at  the  inlet  face.  This  is  accomplished  by  running  SWINT  to  the 
inlet  face  and  saving  TAPE17  which  contains  the  final  flow  field 
information.  The  procedure  for  doing  this  with  the  extended  version  of  SWINT 
is  Identical  to  that  for  the  original  version  and  is  described  in  Reference 
2.  The  output  differs  from  that  of  the  original  SWINT  in  two  respects:  (1)  a 
PTO/PTINF  column  has  been  added  to  the  flow  field  output  which  is  the  ratio  of 
the  local  stagnation  pressure  to  the  freestream  stagnation  pressure;  (2)  The 
program  stops  exactly  at  z  ■  ZEND  rather  that  at  the  first  step  greater  than 
ZEND. 


Interface  Program  COWLI 

The  COWLI  program  rezones  the  flow  field  to  lie  between  the  inner  body 
and  the  cowl.  In  addition,  the  inlet  plane  parameters  described  in  Section 
2.2  are  calculated.  This  program  is  applicable  even  when  the  Inlet  or 
portions  of  it  lie  outside  of  the  flow  field  generated  by  SWINT  at  the  inlet 
face.  At  points  outside  of  the  SWINT  generated  flow  field,  freestream 
conditions  are  assumed  to  exist.  It  is  also  possible  to  use  COWLI  to  generate 
a  starting  flow  field  for  external  calculations  downstream  of  the  inlet  lip, 
even  in  cases  where  the  bow  shock  lies  completely  within  the  inlet 

To  run  the  COWLI  program  the  flow  field  at  the  inlet  plane  generated  by 
SWINT  must  be  attached  as  TAPE1I  and  relevant  quantities  in  namelist  INPUTS 


must  be  defined: 
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BNEW,BZNEW,BPHNEW  -  Inner  wall  boundary  of  the  computational  domain  b(<|>,z) 

and  Its  derivatives  b  (<|>,z)  and  b  (d>,z)  .  (See  Figures 

Z  (ft 

1.2) 

CNEW,CZNEW,CPHNEW  -  Cowl  or  shock  surface  description  c(<j>,z)  and  its 

derivatives  c  (<f>,z)  and  c  (<J>,z).  (See  Figures  1,2) 

Z  (p 

IBODY  Controls  conditions  prescribed  along  inner  boundary. 

0  -  inner  body  shape  is  not  changed.  BNEW,BZNEW,BPHNEW 
need  not  be  specified. 

1  -  inner  body  shape  is  changed  to  BNEW , BZNEW , BPHNEW, 
which  must  be  specified.  Interpolated  flow  values  at  the 
wall  are  turned  tangent  to  the  surface  using  an  oblique 
shock  or  Prantl-Meyer  expansion. 

2  -  Same  as  IBODY  *  1,  except  wall  properties  are  also 
assigned  between  the  body  and  shock.  The  outer  boundary 
is  assumed  to  be  a  shock  with  cz  calculated  from  the  shock 
or  Mach  angle  occuring  at  the 

wall.  c  ■  b  +  (c  -  b  ).  and  c  *  b  . 

'  z  z  ’  <|>  Ip 

ICOWL  Controls  the  conditions  prescribed  at  the  outer  boundary 

and  the  type  of  surface  (i.e.,  shock  or  cowl), 

0  -  outer  surface  geometry  is  not  changed  and 
CNEW , CZNEW , CPHNEW  need  not  be  specified. 
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1  -  outer  surface  Is  a  cowl  with  description 
CNEW,CZNEW,CPHNEW.  Interpolated  flow  values  at  the  cowl 
are  turned  parallel  to  the  cowl  surface. 

2  -  outer  boundary  is  a  Mach  surface.  The  user  specifies 
CNEW,CPHNEW  and  CZNEW  is  computed. 

DDZ  Distance  from  the  cowl  lip  at  which  calculation  is 

started.  Specified  only  when  IBODY  «  2. 

RCLUST  Controls  radial  distribution  of  mesh  points.  Default  is  a 

uniform  distribution.  For  other  cases  enter  (r-b)/(c-b) 
for  each  radial  plane  starting  at  the  body  and  moving 
towards  the  outer  boundary.  The  same  radial  distribution 
is  presumed  on  each  constant  <j>  plane. 

AREA  Reference  area  used  in  calculating  induced  load 

coefficients.  Default  is  the  inner  body  cross-sectional 
area  at  the  inlet  entrance  plane. 

IPRINT  Controls  the  amount  of  printed  output 

0  -  print  only  inlet  plane  flow  field  parameters. 

(i.e.,  Eqs  (4)  and  (5)) 

1  -  (default)  IPRINT  ■  0  output  plus  final  rezoned 


flow  field 
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2  -  IPRINT  ■  1  output  plus  initial  flow  field  from 
Tape  11. 

3  -  IPRINT  -  2  output  plus  Jump  subroutine  messages. 

To  apply  the  interface  program  to  an  inlet  geometry,  options 

IBODY=*0/ I COWL” 1  or  IB0DY”1/IC0WL”1  are  used  depending  on  whether  the  innerbody 
slope  or  surface  is  discontinuous  at  the  inlet  face  plane.  Three  additional 
modes  of  operation  are  IBODY-2/ICOWL-O.IBODY-1/ICOWL-2,  and  IBODY-I/ICOWL-O. 
These  options  are  designed  to  facilitate  restart  of  an  external  calculation 
downstream  of  the  inlet  lip.  The  first, IB0DY”2/IC0WL“0  is  applicable  when  the 
bow  shock  lies  completely  Inside  the  cowl.  Slightly  downstream  of  the  inlet 
lip,  the  flow  values  obtained  by  turning  the  free-stream  tangent  to  the  COWL 
outer  surface  provides  an  estimate  of  the  local  flow  field.  An  alternative 
approach  for  handling  this  situation  which  is  also  applicable  when  the  bow 
shock  is  only  partially  within  the  inlet  is  accomplished  with 
IB0DY”1/IC0WL”2 .  Here  the  outer  edge  of  the  computational  domain  is  defined 
to  be  a  Hach  surface.  As  the  SWINT  calculation  proceeds  downstream  from  the 
lip,  the  shock  or  Prandtl-Meyer  expansion  generated  at  the  body  surface  by 
IBODY-1  propagates  into  the  flow  field  and  merges  with  the  outer  Mach 
surface.  The  final  option,  IBODYal/lCOWLaiO  duplicates  the  function  of  the 
INLET  subroutine  of  SWINT. 

The  output  from  COWLI  consists  of  Tape3  which  is  the  restart  file  for 
SWINT  and  printed  data.  The  amount  of  flow  field  Information  printed  is 
controlled  by  the  parameter  IPRINT  and  the  output  flow  field  quantities  are 
designated  using  the  same  headings  as  found  in  SWINT.  The  items  printed  under 
the  heading  "inlet  plane  flow  field  parameters"  are  described  in  Section 
2.2.  The  induced  force  coefficients  are  followed  by  a  value  labeled  force 
error.  This  number  is  the  percent  dlscrepency  obtained  by  calculating  the 
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forebody  loads  using  direct  pressure  integration  as  opposed  to  Equations  (5). 


Applying  SWINT  to  Inlets 

The  extended  version  of  SWINT  is  applied  to  an  inlet  configuration  in  a 
aanner  similar  to  that  described  in  Reference  2  for  external  configurations. 
Several  additional  variables  must  be  prescribed  along  with  a  description  of 
the  geometry  of  the  cowl. 

The  cowl  geometry  is  described  in  a  manner  analogous  to  that  used  to 
describe  the  body.  The  quantities  c,c  ,c  ,c  ,c  .and  c  .  must  be  specified 

Z  ^  ZZ  Zip 

using  fortran  statements  Inserted  at  the  indicated  locations  in  subroutine 
COWL  of  SWINT.  As  an  example,  consider  a  circular  cowl  starting  at  z  ■  1  with 
a  radius  of  2  and  an  outwards  angle  of  7'  relative  to  the  missile  axis  for 
z  >  1.  The  necessary  statements  describing  the  cowl  are: 


c  -  .12278456 
z 

c  *  2.  +  (z-l)$c 

% 

c  -c  -c-c-0 
zz  tty  $$  4> 

The  final  statement  is  not  required  since  default  derivative  values  are  0. 

The  additional  variables  which  control  the  computation  of  the  cowl 
surface  are  specified  in  namelist  INPUT1: 


ICOWL 


0  if  the  outer  boundary  is  a  shock  and  1  if  the  outer 
boundary  is  a  wall. 
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ICOWOPT  0;  cowl  slope  and  surface  Is  continuous  at  starting  plane. 

1;  A  cowl  slope  or  surface  discontinuity  occurs  at  the 
starting  plane.  This  results  in  the  cowl  surface  normal 
derivatives  being  modified  for  4  steps. 

MOD1C  Controls  order  of  accuracy  of  cowl  boundary  conditions. 

0  -  first  order,  1  -  second  order. 

ISWSMOC  Controls  application  of  entropy  extrapolation  at  cowl. 

Extrapolation  is  applied  on  planes  M< ISWSMOC 

ISWMODC  Controls  form  of  cowl  boundary  equations: 

0  -  Form  2A  and  3A  (Analogous  to  14A  and  15A  of  Ref.  2  for 
centerbody) 

3  -  Form  2C  and  3C  (Analogous  to  14C  and  15C  of  Ref.  2  for 
centerbody) 

The  variables  M0D1C,  ISWMODC  and  ISWSMOC  are  analogous  to  the  body 
variables  MODI,  ISWMOD  and  ISWSMO  respectively.  Recommended  values  for  these 
parameters  are  discussed  in  Reference  2. 

The  output  from  a  run  for  an  inlet  configuration  differs  from  an  external 
flow  field  calculation  in  the  following  respects: 

1.  At  each  plane  where  the  flow  field  is  printed,  PTAVR/PTINF  is 
calculated  which  is  the  area  weighted  average  of  the  recovery 
pressure  divided  by  the  free-stream  recovery  pressure. 
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2.  At  the  completion  of  the  calculation,  the  cowl  pressures  are 
printed. 

3.  The  force  and  moment  coefficients  represent  the  integration  of 
pressure  over  both  the  centerbody  and  cowl  surfaces,  but  do  not 


include  induced  loads 


APPENDIX  D.  SAMPLE  RUN 


This  appendix  illustrates  the  application  of  the  extended  SWXNT  code  to  the 
inlet  configuration  shown  in  Figure  D-l.  The  initial  flow  field  was  calculated  at 
z  =  .1  using  the  approximate  conical  starting  program  START  which  is  described  in 


Reference  2.  The  data  cards  used  for  this  run  are: 

SINBUTS 
NC  «  19. 

NO  «  13. 

AON  -  3.3, 

ATTACK  ■  3.0, 

Bill  ■  0.017633, 

Z  »  0.1. 

ZS  ■  0.1. 

SENO 


This  program  generates  TAPE3  which  is  the  SWINT  starting  tape. 

Using  the  extended  version  of  SWINT  with  the  errata  update  file,  the  flow  field 
was  marched  from  the  starting  plane  at  z=.l  to  the  start  of  the  cowl  which  is  locat¬ 
ed  at  z  =  3.216.  The  errata  update  file  is  as  follows: 


•IDENT  ERRATA 
■DELETE  E06E.134.EDGE.135 
DUN  »  0.0 

CALL  JUHPF  (II.  N,  M.  DUN1 

IF  (MN  .NE.  Ml  CALL  JUMPF  (12,  N,  MM.  DUM) 

•DELETE  EDGE. 7S 

1  ABS  (PHI (Nl  -  PHKMPII  •  0.3S  /  PINFI 
•INSERT  FRINGE. 18 

8(N1>  >  B(M2) 

BZ(N1)  •  BZ ( M2 I 
BZZ (Ml I  •  BZZ (M2) 

BZPHI  (Ml)  •  BZPMKM2) 

BPHl(Ml)  *  BPHKM2) 

BPHPHKM1I  •  BPMPMI  (M2) 

•DELETE  REZ0NE.4S. REZONE. 46 
•INSERT  REZONE. 40 

OX  •  1.0  /  FLOAT  (NA) 

DDX  •  NA 
•INSERT  FIEL0.3! 

IF  (IFIN  .CO.  0)  60  TO  16 
•INSERT  FIELD. 36 
16  CONTINUE 

•OELETE  INTES.41.INTE6.43 

DTHETA  ■  TH ( 12* 1 )  -  TM(I.l) 

IF  (M  •  (1  -  IDYAU)  .£0.  21  0THETA  •  2.0  •  TM(I.l) 

IF  (M  •  (1  -  IDYAH)  .E0.  MP1)  DTHETA  ■  2.0  •  (PHIO  -  THII.ll) 
•DELETE  TRAN6D.27.TRAN80.30 
B  CONTINUE 

DY  *  1.0  /  FLOAT  (NSGD) 

S6DM2  •  SGDINSeO)  •  1.0 
SSO(ll  ■  S60(NS60P1>  -  1.0 
SGOP2  *  1.0  *  SSO ( 3 ) 

S6DINS6DP2)  ■  1.0  *  560(2) 

•OELETE  SHOCK. 68 

DCUZ(I)  ■  (FAC1  •  OOX  «  FAC2  •  DDY  .  CES(I>>  /  PINF 
•OELETE  FRINGE. 35 

PYY  ■  (P(N.N)  -  PIN.MPII  /  (PMI(M)  -  PHI  IMP) )  /  PINF 
•INSERT  FRIN6E.26 

R(N.M1)  ■  RIN.M2) 
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•DELETE  OUT. 208. OUT. 211 
CYaXK0»FY 

CHXsXK 1  * ( MX  *  ZC»FYI 
CNY*XK1* (NY  -  2C«FNI 
CMZ«XK1*M2 

•DELETE  OUT. 214, OUT. 215 

IF (CN.NE.01  XCPP«ZC/ZREF  «  CHY/CN  ♦  ZO/ZREF 
IF (CY.NE.01  XCPY»ZC/ZREF  -  CMX/CY  ♦  ZO/ZREF 
•DELETE  FEVAL.41 

SF (61 *RF  • (COSPP*SF (SI »S INPP»SF (111 
•DELETE  REZONE. 60 
•OELETE  REZONE. 82 


The  configuration  geometry  was  described  using  the  following  update  deck 


center-body: 

•IDENT  PRESSLY 
•S  BODY. 19 

REAL  XHCUOOI  t  YHC  (1001  .ZSAVEOI  ,BSAVE(3I 

DATA  XHC/0.0,  4.0.  4.1.  4.2,  4.3.  4.4,  4.5,  4.55.  4.6.  4.65. 

1  4.7.  4.8.  4.9,  5.1.  5.3,  5.5*  5.6,  5.7,  5.8,  5.9,  6./ 

DATA  YHC/0. ,  .70532.  .7228,  .7387.  .7512,  .759.  .7625, 

1  .763,  .7625,  .7811,  .7585,  .7504,  .7391,  .712,  .6829,  .6525, 

2  .6362,  .618,  .5973,  .5744,  .5467/ 

OATA  KK  /0/ 

NPTS-20 
KMKK.l 
ZOIFFa.001 
IF(KK.8T.2IZDIFF«0Z 
ZSAVE (II »Z-ZDIFF 
ZSAVE (21 *2 
ZSAVEOI  aZ*ZDIFF 
DO  30  J-1.3 

DO  20  I-1.NP7S 

IF ( (XMC (II-ZSAVE(JI1* (XHC ( I*11-ZSAVE(J1I.GT.O.)GO  TO  20 
FAC* (ZSAVE ( J I -XHC( 1 1 1 / (XHC ( I • 1 1 -XMC( III 
BSAVE ( Jl «YNC ( 1 1  * (YHC ( I* 1 1  -YHC 1 1 1 1  »FAC 
80  TO  30 
28  CONTINUE 

HRITE (6.20001 Z, 20 IFF, ZSAVE 

2006  FORNAT (•1GEOMETRY  OUT  OF  RANGE  Z.ZOIFF, ZSAVE  •.5E16.B1 

STOPMBAD  OEOHtTRY" 

36  CONTINUE 

B8Z«(BSAVE(3I-BSAVE(1II/(2.«Z0IFF1 
BBZZ»(BSAVE(31-2.»BSAVE(2I,BSAVE(11 I/(ZDIFF«ZDIFF1 
•I  BODY. 35 

IF (Z.LE.4.160  TO  18 
B(NI-BSAVE(2I 
B2 INJ«BB2 
BZZ (Ml aBBZZ 
80  TO  19 

18  CONTINUE 

B2 (HI ». 17633 
B(H|sZ«8ZIHI 

19  CONTINUE 


cowl : 

•B  C0VLADD.56 

REAL  XHCUOOI  .YHCIIOOI  .ZSAVEOI  .CSAVEOl 

DATA  XHC/2.656,  3.1,  3.2,  3.4,  3.6,  3.8,  4.0,  4.1,  4.2,  4.25, 

1  4.3,  4.4,  4.5,  4.55,  4.6,  4.65,  4.7,  4.8,  4,9,  5»,  5.1,  5.6* 

2  5.8,  5.9,  6.0/ 

DATA  YMC/1 . ,  1.004188,  1.0054,  1.0051,  .99996,  .9882,  .9681, 

1  .954,  .9364,  .9261,  .9154,  .8949,  .8768,  .8695,  .864, 

2  .86,  .8572,  .8533,  .8511,  .8502,  .85,  .85,  .6574, 

3  .8646,  .8735/ 

OATA  KK/0/ 

NPTS  *24 

KK-KK.l 

ZOIFFa.OOl 

IF (KK.BT .2) ZDIFFaOZ 

2Z-Z-.356 

ZSAVE ( 1 I aZZ-ZDIFF 

ZSAVE (21 aZZ 

ZSAVE (31 aZZ*ZD!FF 


DO  30  J*1.3 

DO  20  1*1 .NPTS 

IFI<XMCm-ZSAVE(J))*(XMC(l«ll-ZSAVE<Jl>.GT.O.)GO  TO  20 
FAC*  IZSAVE  ( J)  -XMC » 1 1 )  /  ( XHC  ( I  ♦  1 1  -XMC  ( 1 ) ) 

CSAVF  < J) =YHC( 1 > • (YMC< I»> )-YHt( II >*FAC 
60  TO  30 
20  CONTINUE 

UNITE <6.20001 Z.ZOIFF.ZSAVE 

2000  FORMAT (•16E0NETRY  OUT  OF  RANGE  Z.ZDIFF.ZSAVE  •.5EI6.0) 

STOP"BAO  COWL  GEOMETRY" 

30  CONTINUE 

CCZ*(CSAVE(3) -CSAVE <111/<2.*ZDIFFI 

CCZZ*(CSAVE(3I-2.*CSAVE(2)»CSAVE(1>)/(ZDIFF*ZD1FF) 

•I  C0DLA00.73 

C«M)«CSAVE(2i 
CZ  INI  *CCZ 
CZZINIaCCZZ 


The  cards  designated  as  COWLADD.73  and  COWLADD.56  are  so  marked  in  Appendix  B. 


The  namelist  inputs  used  to  make  this  run  were: 


SlNBUTl 

KA*2000.ZENO*3.216. 

ISNNODC*O.HOD1C*O.ISMNOO*O.M001*0. 

SENO 

SOUTRD 

KOUT 111*40. 

SEND 


Sample  output  sheets  from  this  run  are  shown  in  Table  D-l.  The  restart  file  from 


this  run  is  written  to  TAPE17. 


To  complete  the  calculation  of  the  inlet  configuration,  program  COWLI  is  used 


to  rezone  the  flow  field  so  that  it  lies  within  the  inlet.  To  run  COWLI,  the  restart 


tape  generated  by  SWINT  at  z  =  3.216  is  accessed  as  TAPE11  and  namelist  quantities 


must  be  defined.  The  data  cards  used  for  this  run  are: 


S INPUTS 

ICO«L*1,CNEN*1.<CZNEX*.017A,CPHNE**0.« 

IPNINT*0, 

SENO 


The  output  from  COWLI  is  shown  in  Table  D-2.  The  rezoned  restart  file  generated 


by  COWLI  is  TAPE3. 


The  inlet  section  of  the  configuration  is  run  using  the  restart  file  generated 


by  COWLI  which  is  accessed  as  TAPE3.  The  data  cards  used  to  complete  this  run  are: 

$ INPUT  1 

ICOMLsl , 1C0W0PT*] , 

KA*2000«ZENO*S,2, 

I SMMOOC*0  .MOO  1  C*0 .  ISXM00-0.MO0 1  *0 , 

SEND 
SOUTRO 
KOUT 111*40. 

SEND 

The  output  from  this  run  is  shown  in  Table  D-3. 
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Figure  D-l.  Inlet  Configuration  Geometry 
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